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. . P. H.E. F A 0„E 

The work here presented was au^^ested to me by Dr. 
J. Ackeret, and was carried out at the Institut fur 
Aerodynamik der E.T.Hi Problems in the field of super- 
sonic flowe occur with Increasing frequency in recent 
times. It is of interest firet to investigate as to how 
far the relation' extends between the flow of a liquid. on 
a horisontal bottom with the two-dimensi onal flow of a 
compressible ?as. Secondly, problems in the field of 
water flows may be solved directly by the methods of the 
theory of ?as dynamics* which, in r.ecetit years, have "been 
highly developed. 

The present work was undertaken with two objects in 
view. In the first place, it is considered as a contri- 
bution to the water analogy of fijas flows, and secondly, a 

lars;e portion is devoted to the general theory of the two- 
dimensional supersonic flowB. An attempt has been made 
to brin*; the latter Into such shape and detail as to facil- 
itate as much as possible its application by the engineer, 

who is lees familiar with the subject. 

Here, I should like to sxpreas my thanks to Dr. 
Ackeret for his encourcqement nnd aid, and to Dr. de Hnller, 
Assistant at the Institut fur Aerodynnmik, for his friendly 
support . 



Translator's note: The term "fijas dynamics" is defined in 

the Introduction. 
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APPLICATIOB OF THE METHODS OP OAS DYNAMICS TO 
WATER FLOWS WITH FREE SUBPACE* 
PART I. FLOWS WITH NO ENEHflY DISSIPATION** 
By Ernst Preiswerk 

IBTRODUCTIOB 



Let there be considered a gas at rest in space or a 
portion of space, and let a piston or a movable portion of 
the boundary set the gas in motion. In the case of an in- 
compressible fluid, the latter will begin to flow simulta- 
neously over the entire space at the Instant the disturb- 
ance is applied. With a compressible fluid the case is 
otherwise. The effect of a disturbance first shows up in 
a restricted portion of the space only at a definite time 
interval after the start of the disturbance. If the latter 
iB small, the speed of propagation of its effect is equal 
to the velocity of sound in the qas . In an ideal gas. it 
is proportional to the square root of the absolute temper- 
ature T and depends only on the latter. 

If the velocity of flow in a fluid is small compared 
to the velocity of sound, the fluid may be treated as in- 
compressible . The relation between velocity c (m/s) 
and pressure p (kg/m a ) at various pofnts of the flow, is 
in the case of absence of friction, given by the Bernoulli 
equation. As soon, however, as the velocity differences 
at various points of the flow attain the order of magni- 
tude of the velocity of sound, the compressibility of the 
gas may no longer be neglected. Density P (mass per 
unit volume, kg s a /m*) and temperature are variable, so 
that the laws of thermodynamics must be taken Into account. 
Thethoory of such flow comes under Gas Dynamics (refer- 
ences 1 and 7) . 

*"Anwendung gasdynamischer Methoden auf Wasserstromungen 
mit freior Oberf lache. " Mitteilungen aug dem Insti- 
tut .fur Aerodynamik, M No. 7, 1938, Eidgenossische Tech- 
nische Hochschule, Zurich, 



**For Part II, see N.A.C.A. Technical Memorandum No. 935 
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Depending on whether the flow velocity is smaller or 
larger than the velocity of sound, we speak of a subsonic 
and a supersonic flow, respectively, the two kinds being 
e ssentially . different in character. They may occur side by 
side in the same flow since the velocity c and the sound 
velocity a in general vary from point to point. The quo- 
tient velocity of flow per velocity of sound for a definite 
point of the flow is denoted as the local Mach number U = 
c/a (roforence 4) . For K < 1 the flow is subsonic: il 
> 1, supersonic. The subsonic flows in the neighborhood nf 
M = 1 have as yet been little investigated. To are far 
bettor acquainted with the properties of supersonic flo-»s, 
though chiefly the two-dimonsional flows!* 

Between the variables, pressure, temperature, and den- 
sity, there holds the equation of state for an ideal gas 

P = g H p T / / 

where H. (kg m/kg° = m/o) is a .constant that is different 
for en.cn gas. By the addition of hoat , compression, and 
expansion, all possible states may bo attained in the gas. 
If, however, hoat is neither added nor taken away, and in 
the ^as itself no heat arises through friction then, in 
addition tn equation (l), the following adiabatic equations 
hold b^twoon the state variables: 

p/p o = (P/P 0 (2a) 

PIP, = (T7T 0 j' l/K - 1 (2b) 

P/P 0 = (T/T^/*- 1 (2c) 

where p 0 , p 0 , p D is any reference state, and k is con- 
stant for an ideal gas, being the ratio of the specific 
heat at constant ijrensure (c v ) to the specific heat at 

constant volume (oy)'. This case of adiabatic change of 
stato is the one thnt obtains in an ideal flow (no fric- 
tion, no addition of heat from the outside, heat conduc- 
tion and heat radiation in the flow itself negligible). 
As reference state in a flow thero is generally chosen the 
state at a point of rest. 

In order to be able to apply readily the energy equa- 
tion to thermal processes., there is introduced a further, 

*1) Three-dimensional flows: references 6. 8.. 20.. 26.. 29. 

2) Two-dimensional flows:, references 1 (or 2), (pp. 308- 

322); 3, 7 (pp. 407-444), 14, 15, 17, 18, 27. . 

3) Transition region of subsonic, and supersonic flows: 

references 9, 14 (pp. 57-67), 28, 30. 
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state variable, namel , the heat content i, defined by 
i =. CpJE .('Injce m/kg 3*. Let the heat content at a point 

of reat be 1 Q . The flow velocity at an arbitrary point 
(i , P, T , p) of the flow Is then computed from the energy 
equation to be 

c B = 2« (i 0 - i) = 2* c p (T 0 - T) 
Transforming with the aid of equations (l) and (2) 

--^['-sTi 

This equation (jiVea the relation "between the pressure and 

velocity for the compressible adiallatic flow and replaces 
the Bernoulli equation. To a first approximation, i.e., 
for small Mach numbers, it *joea over Into the Bernoulli 
equation. Bor the velocity of Round, ma have 

a E = dp/dp (reference 13, p.536)(4) 

or, using equation (2a): 



(3) 
(3a) 



a a - k £ =s;lfR T (4a) 
P 



From (?a) and (4a) there is obtained: 



M = c 



k-i 
k 



Jrom the adiabatic 'equation (2a) 



PO p " J \v 



k = , p°- k 



*The heat content is usuall;' expressed in kcal/kg, Many 
computations are sijnplif ied, however, if the heat is con- 
sistently expressed in mkq instead-of kcal. The specific 
heats cp and Cy must then be ^iven in mk?/k§° instead of 
in kcal/ks°. The carrying alonf; of the faotor A = 1/427 
k?m/kcal is thereby avoided. H is simply Cp - c v , etc. 
In what follows, this assumption will everywhere be used. 
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and substituting in the above equation and solving for p Q 
we have 



PO 



Expanding the brackets into a aeries there is obtained : 
P o = P L 1 + feU — M k^l Vk^l "" / 1x2 \~Z~ ) + * ' ' 

Po - P - P [A M * + —] 
a k 

The common factor M can be taken outside the brackets 



p -p=p S M 3 + i H a + li2-ki 4 l(2-k)(3-2k) M 

Consider 

P 0 » m P °! a a 

2 2 a a 

Substituting a a from equation (4a) : 

| 0° = K a k | 
We thus have, finally 

p o -P--seU+ 3 M + -, gt / M + .... J 



6 + 



(5)1 



For M "0, the above becomes the Bernoulli equatiorf 

c a = po - p. A better ap-oroximati on is Ij c a = (p c - p)/ 

(1 + i M S ). The first two coefficients, 1 and l/4, in the 

series are independent of k. For k= 1.4, the next two 
coefficients are l/40 and l/l600. 

We shall now "brin? out an important property of the 
supersonic flows. Let us consider first a parallel flow 
with constant velocity c. The velocity of sound corre- 
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sponding to the temperature of the gas also hae the same 
value over the entire flow plane. ' If a small cylindrical 
obstacle is situated in such a supersonic flow, the dis- 
turbance produced "by the obstacle is propagated with re- 
spect to the. moving gas with the local sound velocity. 
The waves are circular cylindrical in shape (fig. l). Let 
the obstacle be located at point P. If the wave center 
K x is at point X, a time interval t => z/o, has passed 
since this wave arose. It then has the radius r = a t = 
a x/c. At the point F such waves arise continuously. 
All of them have as their oommon envelope two straight 
rays, the Mach rays, which form with the direction of flow 
the Mach angle a; sin a = r/x = a/o. If the obstacle at 
P is small, the intensity of the circular waves is small 
to a higher order. Only along the liach rays are the circu- 
lar waves dense enough for the effect of the disturbance 
to be of the order of magnitude of the lattor. The effect 
of a disturbance at P is propagated only along the Mach 
rays through P. Now instead of a parallel flow, we shall 
consider a general supersonic flow. The flow velocity and 
the sound velocity vary from point to point. Por each 
sufficiently small partial region of flow tho .same consid- 
erations as above are valid, tho direction and Mach angle 
varying only from point to point. The disturbance arising 
from a small obstacle at P is now propagated along curved 
linos (fig. 2), those being Icnorn as Mach linos. Por each 
flow there are two familios of Mach linos. All effocts 
arising from the boundary of tho flow are evidenced along 
those lines of the flow. 

It is possible with liquid fluidB (water) to produce 
flows that show a far-reaching analogy to the dimensional 
flows of a compressible gas (references 5, 11, 13 (p. 537), 
21, 22, 23, and 24). 

A flow of this kind is obtained if water is allowed 
to flow over a horizontal bottom under the effect of grav- 
ity. The surface of the water is assumed to be free. At 
the sides it must be bounded by vertical walls or it must 
flow into water of a definite depth at rest. The fixed 
vertical walls correspond to the boundaries of the gas 
flow. A channel with horizontal bottom and rectangular 
cross section with variable, width, the axis of which need 
not be rectilinear, is an example of this type of boundary. 
The water ■ flowing into water at rest corresponds to a free 
gas jet. An open sluice, from which the water flows out, 
is an example of the second boundary condition. The bottoms 
of the upstream and downstream water must lie in the same - 
horizontal plane. 
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She velocities that occur in such flows are very 
8 mall in comparison with the sound velocity in water 
(about 1,430 m/s). The latter plays no part at all in 
the considerations that follow. It is- another velocity 
which is analogous to the velocity of sound in a gas. 

In the present work only stationary flows will ho in- 
vestigated. The free upper surface of the water is' then 
a fixed surface in space. The water depth h varies from 
point to .point of the flow. For each depth there exists 
for Ion? plane waves a wave propagation velocity </gh, 
which depends on the depth alone. On the basis of this 
wave velocity the water flows described may be divided into 
two groups which, as in tho r case of the gases, differ es- 
sentially in their properties. If the water velocity is 

less than -s/gh", the water will be said to "stream"; if 
greater than */gh, the rater will bo said to "shoot." 

PAHT I. PLOWS UITH NO ENER&Y DISSIPATION 

Differential Equation of the Water Flow 

1. Energy Equation 

It will be assumed that the flow of the water is fric- 
tionlens so that conversion of energy into heat is exclud- 
ed. The onergy equation then simply states that the sum of 
the potential and kinetic energy of a water particle is 
constant during its notion. 

Let us consider a flow filament (fig. 3) which passes 
through the point y 0 , z 0 of the initial cross section 
x = 0. Along this filament, between the pressure p and 
the velocity c, there obtains the Bernoulli equation 

.p + | c a + p g z = p x "+ | Cj a + p g a,, . (6) 

On tho surface of the water p is constant and equal to 
the atmospheric pressure p-g. In what follows we may, 

without error, set this equal to zoro since only pressure 
differences are of physical significance in tho case of 
incompressible flows. Tho magnitudes denoted with the 
subscript 1 refer to an arbitrary but fixed point of the 
flow filament (reference point); The magnitudes without 
subscript ' refer to a variable point. If the water flows 
out from an infinitely wide basin, then the velocity in 
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the "basin is o Q = 0. Also, the curvature of the free sur- 

face is zero. The plane x = 0 is assumed to lie in this 
region. We choose the point x 0 , y Q , z 0 . as reference 
point. The corresponding water depth is denoted by h Q 
and is at the same time the maximum depth occurring. 

Tor the above reference point, the Bernoulli equation 
reads : 

P a ' A 
p + -j5C+pgss=p 0 + pgz 0 

from which 

c a = 2g(z 0 - a) + 2( Po - p)/p (7) 

Te now make a simplifying assumption, namely, that 
the vertical acceleration of the water is negligible com- 
pared with the acceleration of gravity. Under this assump- 
tion the static pressure at a point of the field of flow 
depends linearly on the vertical distance under the free 
surface at that positibtf: 

P 0 = P «< h o - a o> (8a) 

and 

p = p g(h - z) (8b) 

The above substituted in (7) ^ives, finally, 

c a = 2g(h 0 - h) = 2g Ah (9) 

The energy equation (9) holds for the flow filament pass- 
ing through y Q and z 0 at x = 0. Since, however, at 

x = 0, all the stream filaments that lie one above the 
other, have the same h 0 and for all of them, c 0 = 0; 
and since equation (9). does not contain the coordinate z, 
the velocity c at x, y . is constant over the entire 
depth and is given only by the differenco in height Ah 
between the total head and the free level, Ah being, at 
most, equal to h Q . The maximum attainable velocity there- 
fore is c mai = J 2g h Q . The energy equation' may thus "be 
written 

(c/c nax ) a = c a /2g h Q = Ah/h 0 (9a) 



In a gas the maximum velocity is c Bax = J 2g i 0 , 
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and equation (3), corresponding to (9a), becomes: 

(c/cinai) 8 = c 8 /2S i 0 = Ai/i 0 = AT/T 0 (10) 

From these two equations it may "be seen that the ratio of 
the Telocity to the maximum velocity for the water and gas 
flows "becomes equally large if 

(h 0 - h)/h 0 = (T G - T)/T 0 

This is the case for 

h/h 0 = T/T 0 

With respect to the velocity there exists therefore an anal- 
ogy "between the two flows if the depth ratios h/h 0 are 
compared with the gas-temperature ratios T/T Q . The water 
depth corresponds to the gas temperature, and conversely.* 



2. Equation of Continuity (reference 1?, p. 320) 

TTe shall set up the equation of continuity in differ- 
ential form. For this purpose we consider at x, y a 
small fluid prism of edges dx and dy and height h (fig. 
4). Let u and v he the horizontal components, and w 
the vertical component of the velocity c in the direction 
of the coordinate axes x, y, and z. 

Neglecting the vertical acceleration of the water in 
comparison with the acceleration of gravity, equation (8b) 
is valid. From it, we have: 

dx r ^ 3x By 3y 



The right sides of the above relations are Independent of 
z, so that the horizontal accelerations for all points 
along a vertical also are independent of z. The horizon- 
tal velocity components u and v are thus constant over 
the entfLre depth "because they were so in the initial state 
(of rest). 



It is not a question of setting absolute values of the 
velocities equal to each other "but only, of course, non- 
dimensional magnitudes, as c / c max a 
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The continuity- equation for the stationary flow sim- 
ply expre a a-es the f act - that the- quantity . of, fluid flowing 
into the prism (fig. 4) per unit time is equal to the out- 
flowing mass. Since the density pf the water is constant, 
the same holds true for the inflowing fluid volume dq e 
(m 3 /s) and for the outflowing volume dq a ; dq e = ia . a * ' In 
the x-direotion the volume u h dy enters per unit time; 
dq Q "becomes = u h dy + v h dx. The total outflowing vol- 
ume, except for infinitely small magnitudes of higher order 
"becomes : 




This continuity condition written out and divided by dx dy 
gives the equation of continuity 

B (h u) + 9( h v) _ 0 ' ' (11) 
8x By. 

The continuity equation for a two-dimonsi onal compressible 
gas flow is 

u > + 5(P v) = o . ( 12 ) 
9x 3y 

Comparison of the two equations (ll) and (12) shows 
that, just as the energy equations., the equations of con- 
tinuity for the two flows have the same form. From these 
we may derive a further condition for the analogy, that 
the specific mass p of the gas flow corresponds to the 
water depth h. It may be clearly seen now why the incom- 
pressible flow of the water may hear a relationship to the 
flow of a compressible gas. As a. consequence of the com- 
pressibility in a two-dimensional gas flow, the gaB mass 
per unit of bottom area is not a constant but varies from 
point to point of the flow plane. Since the water depth 
in the flow, with free surface varies, the mass per unit 
bottom 'area for this flow is also a variable. 

From the equation of continuity t we derived the result 
that the water depth h corresponds to the specific mass 
p. By' comparison of the energy equations of the two flows, 
it followed, however, that the water depth h was simul- 
taneously also the analogous magnitude for the temperature 
T. This is possible without contradiction only if a very 
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definite as sumption is also made as regards the nature of 
the comparison gas. For the gas flow p depends upon T, 
the relation "between the two "being the adiabatic equation, 
(2D) 

P/P 0 - (T/To) 17 ^ 1 

Now p/po = k/h 0 and simultaneously T/T 0 = h/h 0 , and 
substituting in (2b), we* have the equation: 

h/h 0 = (h/ho) 1 ^" 1 

which obviously is satisfied only for 

k = 2 (13) 

Thus we have the result that the flow of the water is com- 
parable with the flow of a gas having a ratio k = Cp/ c v = 

2. Such gases are not found in nature. There are, how- 
ever, many phenomena which do not depend strongly on the 
value of k, so that the analogy, has significance also 
for actual gases. 



3. I rrotat ional Motion 



Before introducing the condition of absence of vor- 
ticity, we make a slight transformation of the continuity 
equation (ll), taking account of the energy equation (9). 
The latter solved for h, reads: 

h = h 0 - c a /2g 

Hence 

8h _ _ _1_ 3(c a ) 
3x 2g 3x 

and using the fact that c a = u a + V s ,* this gives 

|li = - I ( u §S + r (a-) 



Since u and v are constant on a vertical, and since 



from (9), c also is constant, w = J c a - (u a + v s ) is 

also constant, and since w vanishes at the bottom, it 
may be neglected in comparison with the components u and 
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Similarly, 

: ^V"-17 u > + v igN ' Cb) 

ay « \ or ay/ 

The continuity equation (11 ) may also "be written in the 
form 

f3h + !£u + !3h + f^T = o 

9x 3x 3y 6y 

Substituting in the above the ezpreeeions (a) and (b), 
there is obtained: 



3x g \ 3x Bx/ 3y g \ 3y 3y/ 
The above equation divided "by h and rearranged, gives: 

TTe now introduce the condition for absence of vortic- 
ity. This will he true if |Z - a o. In this case, 

ax ay 

there exists a function $(x,y), the velocity potentiali 
of the coordinates x, y such that 

u = a* v = aj» 

3x 8y 

Substituting $(x, y) into equation (14), the latter may 
be written:* 

This is the differential equation for the velocity poten- 
tial of the ideal free surface water flow over a horizontal 
bottom. The equation is partial of the second order and 

— — — — — ■ ■ — — j 

* 

Instead of we write in what follows in the usual 

3x 

notation 3-* = == etc. 

x bx 3x3y 17 
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linear in the second derivatives. The .coefficients depend 
on the derivatives of the first order and on these only. 
It is to '"be observed that g h 1b not a constant but, ac- 
cording to the energy equation is 



gh = gh 0 - c a /2 = g h 0 - ~* g -Z- 



The equation corresponding to (15) for the velocity 
potential of a two-dimensional compressible flow la (ref- 
erence 1 (or 2), p. 308. 

o xx (!-!f:) +477 (x.V)- 2 » iy ^ = o („, 

The two equations (15) and (16) beco me I dentical if 
gh/2gh 0 is replaced by a 8 / 2gl 0 . J g h is the basic wave 
velocity in shallow water, and corresponds to the velocity 
a in the gas flow. 



4 . Summary of the Blow Analogy 

We shall -yet inquire what magnitude in the water flow 

is analogous to the gas pressure. Writing tho equation of 
state (1) for an arbitrary state and for the state at rest, 
there is obtained by division: 

P/P 0 = (P/P 0 ) (T/T o ) 

Substituting for p/p o the corresponding value h/h Q , and 
for T/T q also. h/h 0 , there is obtained the value corre- 
sponding to p/p 0 ! 

p/ Po = (b/h 0 ) 2 (17) 

This is also obtained directly from the adiabatic equation 
(2a) with p/p Q = h/h Q and k = 2. 

The pressure Pq. on the bottom surface is proportion- 
al" to the water depth h; with p^ as specific mass of 

the water pg, = p^fg h. This pressure has no analogy in the 
two-dimensional gas flow. In particular, it is not the 
magnitude corresponding to the ga|| pressure since the cor- 
responding magnitude to p is h and not h. The force 
P of the water flow per unit of length of the vertical 
wall Is, on account of tho linear Increase of the pressure » 
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with distance "below the free surface, given "by 



P = 



Tor P, therefore, we have P/P 0 = (h/h Q ) . Oomparison 
With equation (1.7) showB that p/p 0 = P/P 0 - • Tie: magnitude 
of the water flow' corresponding to the ?as pressure p Is 

thus, the force of the water on a unit strip of the side 
trails.' The pressures in the two-dimensional compressible 
flow are analogous to the forcesln the water on the ver- 
tical wajLlB. 

Prom the differential equation for the velocity Poten- 
tial,, we have derived the fact that the velocity of sound 



Tie dlfferen- 

of the energy 

&*■-*■ ygh is 

consequence of 
of these two 
and for k = 



a corresponds to the wave velocity ,/g h. 

tlal equation arose 'through the comoination 
and continuity equations. Thus the result 
'not something essentially new but Is only a 
the results p -«--»■ h , T ■* — »■ h, and k = 2 
equations. We have a 3 - gkHT = g(k - l)i, 
2 and i-<— -» h , ■ thl s >;I ves a 2 <f vgh. 

Since the velocity corresponding to a is 
there corresponds to the .subsonic flow c/a < 1 the flow 
with c/ ./gh < 1. The water in this case is said to 
"stream," while the water flow corresponding to the super- 
sonic flow is Raid to "shoot." The essential difference 
In character "between the supersonic and subsonic flows ex- 
ists also in the case of water between streaming and shoot- 
ing flows . ■ " 

The analogy considered in this section holds for flows 
with Mach numhers smaller and greater than 1. Essentially, 

however, only the flow of shooting water will be treated in 

this work: Application will therefore be made of the ex- 
tensively developed theory of two-dimensional supersonic 

flows to the flow of water. 
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TABLE OF FLOW ANALOGY 





Two-dimensional gas 
flow 


Liquid flow with free sur- 
face In gravity field 


Natureofthe flow 
jnedium 


Hypothetical gas with 
_ _k = Cp/cv- = 2_ 


Incompressible fluid 
(e.g,. water) 


S{lde boundaries geOj 
metrically similar : 


• ■ "> 


Side boundary vertical 
Bottom horizontal 


Analogous magnltudei 


Velocity c/c max ,c/a* 

Temperature ratio, T/T ( 
Density ratio, p/p 0 
Pressure ratio, p/p Q 

Pressure on the side 

boundary walls 

P/P 0 


Velocity c/Cj^, c/a* 
Water depth ratio, h/ho 
Water depth ratio, h/ho 
Square of rater depth ratio, 
tt/*o> 

Force on the vertical walls. 
P/P 0 = (h/h 0 ) 8 


Sound velocity a 
Mach number c/a 
SuDBonic flow 
Supersonic flow 
Compressive shock 
(right and slant) 


Wave velocity »f gh , 
Mach number c/ ^gh 
Streaming water 
Shooting water 
Hydraulic jump 
(normal and slant) 



MATHEMATICAL BASIS 
5. Introduction. 



For the treatment of fields of flow subjected to the 
boundary conditions, various mathematical methods, depend- 
ing on the type of flow considered, are available. The 
mathematical basis for two-dimensional incompressible flows 
is the conformal transformation method familiar from the 
function theory. For the computation of compressible sub- 
sonic flows, use is made of the theory of general ellip- 
tical differential equations. This theory has not yet been 
sufficiently developed as a practical method. For the com- 
putation of supersonio flows, however, and hence for "shoot- 
ing" water, there haa been perfected the method of charac- 
teristics of the theory of hype r"bo lie partial differential 
equations by Prandtl, Steichen, and Busemann. 



Since the characteristics method is as yet little 
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known and, particularly, since it has not yet been applied 
to the Investigation- of flows of _ ."shooting 11 water, this 
method in what follows, will be di scUBsed "in some detail. 



6. Introduction of New Variables 

The velocity potential $ (x, y) may "^e geomotri oally 

represented by plotting vertically at each point of the 
flow plane x, y the corresponding value of We thus 

obtain a surface in space which we ehall denote as. a $-sur- 
face. The slope of this surface along any direction gives 
the component of the flow velocity in' this direction. 

Let the velocity along a line AS of a shooting flow 
of water be given in magnitude and direction (fig. 5) . 
This velocity at each point of A3 may be decomposed into 
component8 ct and c n , tangential and normal, respec- 
tively, to AB . Simultaneously, there will also be given 
the slopes of the Q-surface corresponding to the flow in 
the two directions and, finally, the value ©(x, y) itself, 
except for a nonessential 'constant, will also be dotermined: 

8 

5 " / T* ds + ®A 

0 

The five magnitudes x, y, $ (point P) and $y 

(slope) are denoted as an element of the Ssurface. An el- 
ement is simply an infinitesimal piece of the ©-surface 
giving the position and elope. The assignment of the ve- 
locity along AB Is equivalent to the assignment of an 
elementary strip of the ©-surface (fiq. 5) . The mathemat- 
ical problem may ■ thus be stated as followe: To find a sur- 
face whose curvature and slope satisfy the differential 
equation ( 15 ) . 

It is possihle to put equation (i5), by a transforma- 
tion of variables, into n. simpler form (reference 27, p.o-10) 

We consider first a ueual coordinate transformation -• 
a so-called "point transformation.™ Let x and y be the 
independent variables; $ a function of x andy,$(x, y). 
Then net variables ■ X, Y, X nag be introduced "by. def ining 
them through the following- equations : 
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X = X(x,y.*(x,y)) 
T = Y(x,y,4&) 
X = X (x,y,$) 



(1-8) 



The function X may "be represented "by a X -surface 
in an X, T, X space, taking X and 7 as ths inde- 
pendent variables. To each point x, y, 4, there corre- 
sponds according to equation (18), an image point X, 7 X. 
Conversely, to each image point corresponds its original- 
point since, in general, equations (18; may he solved for 
z, y, and 4>: 



x = x(X,T, X) 
y = y(X,T, X) 
G = $ (X.T, X) 



;i9) 



Let us, for simplicity, consider first a single inde- 
pendent variahle x and a function 0 = Q(x). The point 
transformation in this case is given "by the two equations: 

X = X(x,$(x)) ' and X=X(x,«) (18a) 

Solving (18a) for. x and $, there 1b obtained: 

x = x(X,X) and $ = &(X,X) (19a) 

To each pair of values x and 4 (point P), there corre- 
sponds according to (18a), a pair of values X and X 
(point P*) (fig. 6). An entire curve has another curve as 
its image and the transformation is uniquely reversihle. 

We shall now consider a more general transformation. 
Let an entire element - that is, x, y, $, $ z , $ y he trans 

formed. In place of formulas (18), we now have the more 
complicated transformation formulas: 



X = X(x,y, fc.^.Sy) 
T = Y(x,y,$,0 x ,$y) 

x = x(x.y.*.$ x -*v) 



(20) 



In the case of a single independent vnriahle, an element is 
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given by the triple x, 4, $ x (p^oint and direction). To 

transform this element the- transof rmation formulas would 
lie' ~ - "* - 

X =* Z(x,*,* x ) and X = X(x,»,* x ) (20a) 

From the above we have : 

dX o Zj. dx + X$ d« + X$ x d$ x ° (X x + X$ * x + $xx) .4* 

and 

d\ = (X, + X$ ® x + X^ fl^x) dx 

4 

go that 

x _ = dX _ *x + Xft $x + *4>i $xx . (21) 

1 dX " X x + X$ » x + Xq x 4 xx 

hence, dx/dX, as(2l) shows, in general depends on x, 

*x«. an< * *xx' If ' for example, 8 curve (fii?. 6) 

is prescribed, then at each point of the curve these four 
values are known. From the three formulas (20a) and (21) 
there are thus determined at each image point P* the val- 
ues X, X, and Xj. There ia thus obtained the curve 
as the intake of curve Correspondingly , may also 

be drawn If the entire curve Xj^ is «5iven. On the other 
hand, from the element x, <t x , It is not possible to 
determine anp element X, X, Xj from the formulas (20a) 
and (21), different elements being obtained, depending on 
how ^xx is chosen. In one case, however, the transfor- 
mation is such that the image of an element is aqain an el- 
ement, and conversely. This is the case when dX/ dX .in 
equation (21) becomes Independent of $ X xi which is true 
only if 

Xx + Xfl * x *$x 

If the transformation formulas (20a) satisfy the condition 
(22), then the elements uniquely correspond to one another 
in the transformation. 

An example of the above is the Le^endre transformation 
of i, $ to X, X', of which, we shall make important use 
"below ; for this transformation, the followinq transforma- 
tion formulas hold: 
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X = » x 

X = ft x x - ft 

We then have : 

dX = ft xx dx 

dX = 0 I dx + ft xx dx x - ^ di = i ft xx dx 

so that 

dX/dX = x, independent of ft xX 

The transformation with corresponding elements has in 
addition, another special property. Let us assume that at 
point' P (fig. 6) two curves ft^ and ft B touch each 
other. They thus have at point P a oommon element x A = 

Xjf ft A = *B- and *xA = *xB : * ut $ xxA =^ $ xxB the clirTeB 

being assumed in contact' but not osculating. According to 
the traneformation formulas (20a) t we shall also have for 
this point. X A = XB and X^ = Xjj. The two imase curves 

Xa and XB then have the point P*. the image of P, 

Also in common. Since,' however, dx/ dX in general, con- 
tains ft IX according to (21), and this second derivative 

is different for the curves A and B, the two linage 
curves will intersect in point P* and not touch as the 
original curves do. Only, if dX/dX is independent of 
$11 will the two Image curves X^ and XB also touch at 

point P*. This is precisely the case for the transforma- 
tion with uniquely reciprocal element correspondence.. For 
this reason such transformations are knownas contact 
traneformatione. • 1 

*1> In correspondence with the concept-point transforma- 
tion, the term "element transformation" is more logical 
than contact transformation. 

2) The transformation (20a) becomes an element transforma- 
tion aO soon as, instead of only the two formulas of (20a), 
three are used: 

x= X(x,*,ft z ) X = XU.ft.ftj.) and X x = Y. z (x r 9,\) ( 2 °t>) 

There then corresponds to each x, ft, ft , an Z ( X, X x , 

and conversely. It is to be noted, ..however, that there is 
a relation between the three variables since Xi = dX/dX. 

If the left' aide of "(20b) is independent of ft XX i the riqht 

side must be. But this is precisely tho contact trans- 
formation. 
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■'The result found above we shall now apply to two inde- 
pendent jaj^aM^s , x, y, and their function ft. The trans- 
formation formula8 are: *" ' 



I = X(x. y. ft, ft x , ft y ) 
T = Y(x, y, ft. » x . ft y ) 
X = X(*. y. ft. ft x . fty) 



(20) 



Since x, 
also ft,. 



Y, and X contain, in addition to x, y, and ft, 
and ft._, there will in" general also occur in 



and 



X 



8X/BI = fj (X,y, ft, 4y , ft^j., fl^y , fcyy) 

ax/ai = f a (x.y.*.« s .4 y .« i:x .4 xy .(e y y) 



(23) 



the second derivatives ft T _, ft-,, , 'b shall interpret 

xx xy yy 

ft(x,y) aO a eurface (fig. 7) . Two surfaces "ft^ and $B i 
which touch at a point, have x, y, ft, ft x , fty in common 

at this point. From the transformation equations they 
will then also have the image point X, Y, X of the con- 
tact point in common. Since, however, Xj nnd. Xy con- 
tain the second derivatives of ft, the tmo transformed 
surfaces will no longer be in contact at the common point: 
(XX), and (XX), not Tieing equal - similarly, (Xy) A 

and (Xy) b " The transformation again gives a unique cor- 
respondence of the elements only If the equations (23) do 



not contain the magnitudes ft. 



ft. 



and ft 



In this 



77" 

go over after 



xx" w xy 

case two surfaces in oontnct at a point, 

tranef ormat Ion Into two surfaces which at the image point 
again have a common tangent plane. 



The Legendre contact transformation for two independent 
variables is 



X =ft 3 
Y =* 3 
X = ft. 



+ fty y - ft 



(24) 



The eurface ft = ft(x,y) with the ahov.e transformation goes 
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over into a surface X = X(X,Y) (fig. -7). We- prove first 
that the above is actually a contact transformation. From 
equation8 (24) 

dX = $xdx + i d$ x + $ydy + y d^ - d4> 

Noting that d$ = $ x dx + 4ydy, three terms drop out. Suo- 
stituting for $ z and $ y , X and Y, respectively, from 
formulas (24 ) , we ha ve 

dX = x dX + y dT 

For the X-surface, the relations must be satisfied: 

dX = XjdX + XydY 

Comparison of the two expressions gives the derivatives of 
X of the first order: 



XX = 1 

Xv = y 



(24a) 



Tnese are independent of the derivatives of $ of the 
second order.. Formulas (24) thus actually express a con- 
tact .transformation, (24) and (24a) /giving the correspond- 
ing element X, T, Xi Xji Xy when the oriqinal element 
i, y, $ x . $y is given. By simple reversal there is 

obtained the element correspondence for the reciprocal 
transformation: 



i = Xx 
7 =xy 

* =XXx+Txy-X 



• x = X 



1 



$y = Y 



} 



(25) 



(25a) 



We wish still to express the derivatives of second or- 
der * xx , $ xy , and * y y 'in the new va r i ab 1 e s X, Y, X. Xj, 

Xyi Xjji Xjy, and Xyy There will then be obtained an im- 
portant result for the applications. 
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Tor this purpose we consider x and y as the inde- 
pendent, variables. Trom the. firsthand second of equations 
(25), there is obtained: 

dZ = XyydX + ^xydY 

dy = Xj^ydX + XyydY 

Solving for dX and dT 

di = ( Xyydz - x XT dy) i/s 
dt = (- x XT dx + x X x<iy) i/n 

where N = (X ZI X T y - ^XY^ 

For the differential of 4, we have ($ -surface) 

dG = $ x dz + flydy (26) 

Substituting in the above (25a), there is obtained: 

di = I dx + T dy 
For the second differential, we have: 

d S 4> = dX dz + dY dy 

a a 

for d x and d y are equal to zero since x and y are 
independent variables. In this equation we. substitute the 
previously found expressions for dX and dT, and obtain: 

d a * = (x T v di a - 2 x XT a* *y + *xx d y a > 

On the other hand, from equation (28) : 

d S 4 = $ xx dx 8 + 2 $ xy dx dy + 4y y dy s 

a a 

Comparison of the coefficients of dz , dy , and dx dy of 
the last two equations shows finally that 

k (27) 



®xx - *YY 1/" 
®yy = *XX l/» 
*xy =-*XY 1/ H 
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These are the required expressions for the derivatives of 
4 of the second order.. 



The coefficients of the differential equation of the 
flow (15) depend on the derivatives of the velocity poten- 
tial 4 of the first order. Introducing new variables 
into that equation (according to the Legendre contact trans- 
formation, the coefficients according to (24) will depend 
on the new independent variables and only on these. The 
partial derivatives of second order will be replaced, ac- 
cording to equations (27), by the partial derivatives of 
second order of the new function with the common denomina- 
tor IT. Since the differential equation (15) is linear 
homogeneous N may be multiplied out. 3y means of the 
Le?endre contact equation, therefore, (15) becomes linear, 
homogeneous, of second ardor, and with coefficients that 
depend on the new independent variables only. 

Let us introduce the new variables X, Y. Physically, 
X and T are the velocity components U and V. The 
new variables according to (24) are: 

(I =0 u = 4 X 

(T =) v = 4>y [ (28) 

X=4 x x+$y7-$=ux+vy-© 

The transformation formulas (25), (25a), and (27) are: 
■ * = X u , y = X v , G=ux+vy-X 

*xx = X vv i/N, * xy = - X^ 1/N, * yy = 1/JT (SO) 

The differential equation (15) in the new variables then 
becomes : 



(29) 



(31) 



x and y beln£ the coordinates of the flow. TTith the . 
Le^endre transformation of equation (15) into (31) , we 
passed from the flow over into its "velocity image" - that 
is, the hodograph (velocity plane) of the flow. At the 
same time, in place of the velo city potential which is 
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a function of the position in the flow, we have introduced 
the "position .determining" potential X,. which is a func- 
tion of the velocity in the hodbgraph. ' "" 

The assignment of the velocity along a curve AB is 
equivalent to the assignment of an elementary strip of the 
©-surface. Since the oontaot transformation is an element 
correspondence, the X-surface must contain the correspond- 
ing X-elementary strip. 

Bor later use, we ehall introduce in equation (31 ) in 
plaoe of the rectangular coordinates u, T, X the cylin- 
drical coordinates c, Cp, X (point transformation), fig- 
ure 8 . 

The new variables are: 




cp = (tan 



(v/u> 



X = X 



whence 



u = c cos cp 



(a) 



= c sin cp 



00 



and 



3c 
3u 



2 



1 



1 



+ v 



a 



2u = cos cp 



3c 
3v 



= sin cp 




3JP _ c o s cp 
3 v c 



We have : 

X = XCu.t) 

so that 



= XCc.cp] = x[o(».t). 3 
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ax ax 3c ax aq> ax m . ax 8 in o 

= + — - s — OO 8 Cp - - t 

Bu 3"c da 39 Bn 3c 3<p c 



(A) 



ax 

a 



i = ax 9c + ax = ax flln o + ax ° 08 9 
t a c av av av a*s °v o 



Furthermore : 



3^X _ 3(3X/3u ) = 3I3X/3ti) c<Jfl _ 3(3X/3u) Bin <p 



3u 



3u 



3c 



jfx. = a(3x/3v) = 3(aX/3v) coa ^ , 3( 3X/3v) sin <p 



<£x = a( 3X/3v) = aOX/av) Bin q, + 3 OX/ 3 v) coa ? 



av 



av 



ao 



acp 



Substituting in the above the values of ax/au and ax/a V 
from equations (A) there is obtained: 

*uu = hoc ™* <P - % ^» + Xcp -^j coa cp 

f« sin cp M cos cp 1 ain cp _ 

- coso? -Xc sincp -Xcpcp— ~* - Hp — ^ J ~ c " 



v a« v 2 sin cp co scp , v,™ sin t> . v sirPro . 
= X cc cos^X^ y + x W H-- + X c 5^- + 



+ Xcp 



2 sine? co scp 




(c) 



and the other two formulas ^ive : 

,2 



v . _ _. v coB a Q3 - sin a cp Y sine? coa cp 
X U v =x cc slnC P coacp+Xcq, * - V? 1 

- x c a lncp cosc ? - Xcp C08a< ?^slnfs ? (d) 

Xvv-Xec ail ^9 + Xoqj * £ + *cpcp £5 *c c 

- xcp 2 ^^^? (e) 
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The transformation formulae (a) to (e) can now be intro- 
duced into i equation (31)-.' The latter then reads in polar 

coordinates: . . ■ , . 

7. Characteristics of the Differential Equation 
(references 10, p. 153, and 31, p. 282) 

The differential equationa (31) and (I51a) are a' spe- 
cial case of the following general form: 

A(X,l) Z XI + 2B(X,Y) Z IY + C(X,Y) Zyy = 

= DiCX.T) Z x + S X (X,T) Zy + ^(X.T) Z (32) 

if for the moment we write Z in place of X. and X and 
Y for u and v , or c and 9, respectively. The coeffi- 
cients A to I 1 of differential equations (32) depend on 
the free variables only. For each pair of variables — 
i.e., for each point of the holograph these three magni- 
tudes are given numbers. There is a simple integration 
method for equation (32) that depends on finding a Taylor 
series for the solution Z = Z(X,Y). 

We seek a solution of (32) that contains a prescribed 

elementary strip. Let the curve, over which the Z-element 

strip is given be expressed in parametric form with t as 
parameter 

X B,Kt^ ) 

(curve A3) 




The Z-surface strip '(the boundary values of Z) over this 
curve is then given by 

Z = F(t) (33) 

and 9Z/3n = Q(t) where n is the normal of the curve AD. 
Along AB: 

dz _ az dx az dT _ az Ylf , az . 
dt - ax It + 5r it - 5x x (t) + 3Y T (t) 
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On the other hand, on account of the prescribed boundary 
"values along the curve AB.we have: 



bo that 



az 

3X ~ " 3T 



X'(t) + || Y'(t) = r i( t ) 



(33a) 



The normal of 'the curve X(t), T(t) has the direction 
.cosines 

cos(n.X) = -Y'(t)/ fvTt) + Y' B (t) 
cos(n.T) = X» (t )//x' B + T» a 



Hence 



|| = ||cos(n.X) + |icos(n.T) = 



dn dX 



3Y 



/x»* -hi 1 ' 



V ax 3y J 



This expression must be equated to G(t). Thus alon^AB 
we also have: 



92 Y'(t) + || X'(t) = J~j7 + G(t) 



ax 



(33b) 



Equations (33a) and (33b) may be solved for 3Z/3X and 
8Z/3Y , since the denominator determinant of the pair of 
equations is 



X' Y« 
- Y ' X 1 

Let the solution be 

az/ax = p(t) 

3Z/3Y = q.(t) 



= X» 8 + Y' 8 ± 0 



(34) 



Differentiating each of these equations with respect to t , 

there is obtained: 
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Z xx I'(t)+ Z XT T«(t)= p ' ( t ) 
Zjy X'{t) + 'Zyy. T« l'('t) 



(35a) 
(35h) 



Tor the second derivatives of Z , we have as third con- 
dition the 'differential equation itself: 



A Zjj + 2B Zjy + 0 Zyy = Di Z x +.BJ. Zy + S*i Z 



(35o) 



If the denominator determinant of the system of equations 
(35) 



X' T 1 0 
0 X 1 Y 1 

A 2B C 



= 0 X 1 8 .- 2B X'Y 1 + A Y 1 8 



(36) 



is not equal to zero, the three equations (35a-c) may be 

solved for Zj X , Zjy, and Zyy. Let there "bo obtained for 

the-derivat ives of Z of the socond order nlongABthe 
values : 



'XX 



= R(t); Z^y = S(t); Zyy = T(t) 



(37) 



Differentiating (35a) and (35"b) with respect to t and 
equation (35c) partially with respect to X and T and 
substituting In the last two equations the values for Z , 
Zj. .. from equations (3?), (34)ar.d (37), there is ob- 
tained the system of equations: 

Z XXX ** 3 + 2Z XXY + Z XYY 7,9 = P" (t) 

Z XXY X ' 2 + 2Z XYY X ' Y ' + Z YYY Y ' 3 = 

= .a(t) 



2 XXX + 2B Z XXY + 0 Z XYY 



A Z XXY + 2B Z XYY + 0 Z YYY 



= 0(t) 



From these equations are obtained the four derivatives of 
third order, of Z along the projection curve of the ?iven 

elementary . strip , since the determinant of the denominator 

is equal to the square of the determinant (36) and thus 

not equal to zero if that determinant is different from zero* 



Proceeding ' in this manner there are obtained all of the 
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higher derivatives of Z starting from the boundary values 
B(t) and Q(t) {equations (33). (34),, (37), etc.). It is 
thus possible to write, the solution of Z = Z(X,Y) also for 
points which do not lie on the curve AB as a Taylor series: 



z(x,Y)=z(xo.T 0 ) + -n [z 1 (x 0 .y 0 )(x-:xo) + M x o.*o>(*- t o> 

+ Yl [Zxx^O.'ToKZ-Zo) 3 + 2Z XT (X 0 .Y 0 )(X-X 0 )(Y-Y 0 ) + 
+ Z TY (X 0 ,Y 0 )(Y-Y 0 ) £ 



+ 

I 



+ • 



This method of solution falls, however, if the determinant 
(36) assumes the value, zero, i.e., if 

o(i,D (H; -aB(x,T)ff ff+ A(I ' T) -o 

or 

C dX S - 2B .dX dY + A dY 3 = 0 . (38) 

This equation, decomposed into linear factors, becomes: 

|A dY-(B+ ( /B 5 -A C) dx| |A dY- (B-.^/b^-A C) dxj= 0 

The denominator determinant (36) thus vanishes if either 



A(X,Y) dY-(B(X,Y)+yB a (X,Y)-A(X,Y) C(X,Y))dX=0 (38a) 

or _ 

A dY - (B - J1F~^ - A C) dX = 0 (38b) 

I't is important tn observe that the pair of equations 
(38a) and (38b) are veil qy the coefficients of the dif- 
ferential equation (32) alone. They are two ordinary dif- 
ferential equations. The solution of each represents a 
family of curves f(X,Y) = k. These two families of curves 
are denoted as the characteristics of differential equation 
(32). If these families of curves, defined by (38a) and 
(38b) are real, then (32) in this region is denotod as hy- 
perbolic. If the two families coincide, then (32) is par- 
abolic. In regions within which the two sets of charac- 
teristics are imaginary, (32) is denoted as an elliptic 
differential equation. 



If, therefore, the curve AB along which the Z-elemen- 
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tarp strip is prescribed as boundary value to (32) is a 
characteristic ,. ..the,., described method of solution by devel- 
opment of Z(X,Y) into a Taylor series fails. 

Ab an application we shall -now compute the character- 
istics of the differential equation of the flow. The com- 
putation is simplest if ve etart from the equation in po- 
lar coordinates (31a) . Comparison of (31a) with (32) shows 
that for this case the magnitudes A, B, and 0 assume the 
following values: 

A.l. B = 0, 0 = - 0 

and the variables X and T' are now c and qp. The ordinary 
differential equations of the characteristics (38a) and 
(38b) then become: 

d<p + y~-±g - l) dc = 0 (39a. b) 

Substituting in the above the energy equation (9): 

*,h = sh 0 - c a /2 

there is obtained the differential equationa of the two 
families of characteristics: 




(40a, b) 



Before we integrate this equation, we wi sh yet to Introduce 
another concept. 

The critical velocity a* (m/e> is ^iven by the con- 
dition that the flow v elo city is -equal to the wave propa- 
gation velocity a =yih, so that the Each number M = 1. 

Thus if c a = z h, a* = c = J % h. Let ub compute the wa- 
ter depth at the critical positions. From the energy 
equation 

c a = 2«h Q - 2gh 
and this should be equal to 
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a 8 = gh 

that is, 



2g h 0 - 2g h = g h bo that h* = | h 0 (41) 



and hence, 



.* _ 



= a 



1 =§«*o < 42 > 



The critical positions in a water flow without energy dis- 
sipation are located where the water depth is two-thirds of 
the total head. These poeitions in an accelerated flow 
are the transition points from 11 streaming" to "shooting" 
water and conversely, for decelerated flow. 

Substituting (42) into equations (40), the latter af- 
ter a small transformation, become: 



1_ /jo/mrf - 1 a(e/a ,) 
/a*) J 1 - (c/a*) B /3 



± dtp = 

(c 

We shall denote c/a* as the velocity ratio c, for which 

a* is taken as the refcronco volocity. Hence, 



1 / c 3 - 1 

° J i - iPTs 



± «P = = / " a d0 (43a, h) 

r 3 

The variables in the above equation are already sepa- 
rated, and the equation map be integrated by a simple 
quadrature. We first introduce a new integration variable i 

— 2 

z = c 

so that we have : 



y (z - i) (3 - z ) 

*T% dz 



y~^~3 + 4z - x S 
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This integral splits up into two parts, J x and J a , of 
which' the first may he directly '.eTOLuated : 



j x = f _ JJ r _Jf__ = /s (sin-) (z-2) 

^ -/ -3+4z-z J J l-(a-2) a 

In -the second integral 

/3 dz 



a - - / a 

Zr-3+4z-z 



we make the anbBtitiition , w = l/ z , so thnt : 

z = l/w 



dz = - -ij dw 



We' now have 



J g = + [* = r___jlggj = ( S in- 1 ) (3w - 2) 

^ 7^3w ff +4w-l ^ 7 l-(Sw-2^ 



= (sin" 1 ) (3/z-2) 
Denoting 



f(c) = / .° Ifl 3 ;. dc (44a) 



1 -l^/S 



we have finally With J x and J a 



f(c) = £ 



J* (sin- 1 ) (c S -2) + (ain- 1 ) (3/ 0 S -2) 1 (44h) 
L J 



The solutions of (43) are thus: 

cp - Cp 1 = f(c) (45a) 
-cp + cp a = f (C) (45b) 
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where <Pi and <p a ■ are the constants of integration -,' 

these being the: parameters of the two families of charac- 
teristics. The latter are shown in figure 9; they are ep- 
icycloids, the loci of the points of the circumference of 

a circle which rolle on another circle (fig. 10). This 
statement can be confirmed in the following manner. 

From the equations (39) (characteristics) . and from 
the energy equation (9), It follows that for h = Othe 
magnitude of the velocity be comes a maximum. in the veloc- 
ity diagram the extremity of c~ a3C then lies on a circle 
Kmax (fig. 9). For all possible velocities that occur, 

c max h > °* For cS > the radicand of (3 9) 

then becomes positive and the root real. Hence, for that 

region of the hodoqraph in which c max > c > «/gh (region 

II), there are two real families of characteristics. This 
holds for the shooting water (supersonic flow) . For a flow 

in which c < ,/gh , the root in (39) becomes imaginary and 

there exist in this region (i)'no real characteristics. 

This is the case for streaming water. 

Let the angle \|f be chosen as parameter (fig. 10) . 
Then, on account of the "rolling condition, " 

a = (r/H) \|/ 

From the triangle PSO , there is obtained for 0 

■ 

0 = (tan- 1 ) \—I »l»±- ] 
{_ (E+r ) - r cos \|/J 

From these two equations, we have: 

cp = a - 0 = (r/E) \|/ - (tan" 1 ) f.^_^sin_jU J ( a ) 

|_(B+r) - r cos \|/J 

From the cosine law for the triangle PTO: 

c = y(H+r) a + r a - 2(E+r) r cos \|/ (b) ' 

Differentiating (a) and (b), there is obtained: 

dip = [(^r) a +r a -2(R4rj r cob M/ J r/E -(E+r)r cos |t r a d>J/ ( c ) 
. (E+r) a + r a - 2(E+r)r cos v|/ 
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do = r < H + r) aln j — dvl/ (d) 

7 "(R+r)" + r* - 2(R+r)r cob \Jf - 

Eliminating in these two equations Bin \|r and cob i|/ 
with the aid of equation (h) , and than dividing (o) "by (d) , 
there ia ohtained:- 

dg> _ 1 e a (a+2r)/H -R(R+2r) (q) . 

do c / 7 - - - -^7^+4^+4^) , H 8 (H+2r) B 

Dividing numerator and denominator of this fraction by 
</c a - R a (H+2r)/H, we have, finally: 



i2 = I / c* - R a ( f ) 

do o J R a - [R/R+2r)] B c 8 

8s was to he proved. Bor R = 1 and (R+2r)/R = Jz , this 
iB the differential equation (43). The epicycloid drawn 
in fi?ure 10 1b thue a charact eri ati c of the family (458). 

The charact eri ati ca of shooting water flow are epicy- 
cloids between two circles whose radii are in the ratio 
</3:l. They are drawn on chart 2 of the supplement. For 
8 gas, the characteristics lie between circles whose radii 

are In the ratio v(k + l)/(k - l) to 1. They are shown 
on chart 1 for air (k = 1.405). 



8. Further Properties of- the Characteristics 

We have seen that If an elementary strip be given as 
"boundary value over the characteristics of a partial dif- 
ferential equation, the solution method by a series devel- 
opment of the required function fails. Gome further prop- 
erties of the oharacterl sties will now be discussed. The 
physical character of the Bupersonlc flow (shooting water) 
which differs essentially from subsonic flow ' (streaming 
aster) - will thereby receive an interesting explanation 
from the mathematical point of view. 

In equation (32): 



A + 2B Zjy + 0 Zyy = Ej_ Zj + E x Zy + Fj, Z 
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let new . variables "be. introduced by making use of a point 
transformation. Let the' new ■ variables 'be: 



A = MX.Y) 
1* = H(X.T) 



(46) 



where for the moment we .-do not fix any definite transfor- 
mation formulas. 'From .(46) we obtain the inverse formulas: 

X = X<> ,n) 

Tho solution of the differential equation (32) Z = Z(X,Y). 

is thus a function of \ and [i. 

Z = Z [X.I*] = Z [>(X,Y),n(X.Y)J 

From tho above, wo havo : 

z y = z x * Y + z^ Mr 



(47a) 



Dif f erentiatlr.^ a second ti:ao, thoro are obtained the de- 
rivatives of second order of Z in the new variables: 

Z XX = Z X*(*X) 8 + 2Z XH X XM-X + Z HH ( ^X) S + Z * *XX + Z H H X X ' 

Z XY = Z XxVl + W^^Vl 5 + WlT + Z * X XY + VxY 

Z TY = Z AX^Y )S + 2Z X^Yl Jt Y + ^n^Y^ + Z ?^YY + Z H**YY 

Putting these expressions in differential equation (32), it 
becomes: 

Z *X [ A *X 8 +2B ^X.V C? Y a ] +2Z Xn [AA x H z +3(X I | iT +X T n I )+GX Tj x Ti + 

+Z^ ^An I a + 2Bn, 2 n T +0ji T a 3 ^aZ^+EjjZ^+ffgZ 

We shall now determine the transformation formulas 
(46). The differential equation of the characteristics is 

0. dX B - 2B dX dY + A dY a = 0 (38) 
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If equation (32) is hyperboli o , (38) has tfo real families 
of curves as solutions.- -Let-. T these be \ _^ 

f, (Z.T) = constant 
and \ ' (49) 

f a (X,T) ="oonstant 

Along each of these curves- . • . ' ' 

f x dX+ fy dY= 0 

This equation together with . (33) si ves for both f x and 

f a , the relation: 

A fj^ + 2B fj fy + C fy" = 0 (50) 

An essential simplification is obtained if, for the 
tranef ormat Ion formulae (46), the following special ones 

are chosen: 



X = f x (X,T) 
p. = f a (X,Y) 



(51) 



[curvilinear coordinatea in the holographs, fig. lib). The 

two coefficients of Z^ and Z^ by (50) then vanish in 

the transformed differential equation, the latter receiving 
the form 



(52) 



This form is called the normal form of the linear hyper—- 
bolic differential equation. It is well suited to numeric 

bal integration by means of the difference method. 

As an application, let the characteristics (45a and b) 

be Introduced as curvilinear coordinates of the position— 
determining potential x (31a). TTe then obtain the nor- 
mal form of the differential equation of flow. 

By elimination of h and h Q from the three equations 

(9) c a = 2$h 0 - 2<sh, (42) a* 8 = 2*jh 0 /3, and a s = «h 

there is obtained: 

c a = 3a* 8 - 2a 8 
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from which, after short computation and substitution of the 
velocity ratio o = c/a*, there is obtained; 

c a 2"g a , , , c a - 1 

T - Ta and S - 1 = 3 " He 

a^-c a 3 - c 

Substi tut in? this expression in (31a) and multiplvina the 
latter by the critical velocity a* (42), then (31a) may- 
be written in nondimensional form: . 

o\ _ o^X 3(c a - 1) _ 5X g (o 3 - 1) _ q 
9c 3 Sep 8 c S (? - c ) dc c(? -» c 8 ) 

In the above we now introduce the coordinates A and H 
through the following expressions : 

x c =*a*c * *\>y>-c 

*cc =x aa (x c) + ^X^cH'c + X^Cprg) + X^cc + *|*»*eo 

3 £ 

X w = Xxx(^ q3 ) + 2Xxn\P^p + X R* (p ^ + *A A cpcp + X H^ 
After substitution and rearrangement , there is obtained: 



a a x 

3- 

oa 



lS>o y c & (5-c ) Wj an" lIjc' c a (3-c a ) \3cp'j 



+ 2 o )t, I dA dp 3(c -1 ) dM- 

dAdM- lcic dc c B (3-c a ) dep 3cp_ 

. 3X jVa _ 3(c a -l) o S X _ 3(0°-!) fl^l . 
. ax Loc c (7-c ' 3cp c(3-c ) o'cJ 

+ ax i"d a u _ ?(c a -i) o a p. _ aCpg-i) a_&] = o (a) 
an LaiF c a (?-^ E ) a^ 3 " c(^-c E ) acj 

The two sets of characteristics (45s) snd (451)) in the im- 
plicit form are now 

f("c) + cp = constant 

f("c) - cp = constant 
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Substituting in (A) for X and p by (51) . the two values 

f 

V-""f(o") +'cp " - - (53a) 

and 

p. = f(c) - qp (53b) 

the coef f icients ' of and X^p, become aero and, since' 

\p = 1, . Pep =~- 1, Ajpcp = 0, Hqxp = 0 

X-5 = df (c)/do = df (c)/dc- ■ 

Xec = d a f (c)/dc a = A a f(*)/ac° 

(A) "becomes: 



8 B 

2 ax ■ 



!Adc/ c a (3-c" S TJ Lax ajij L dc a . ~c(3-c a ) dc x 



= o 



axdn 

and the normal form finally reads; 

d a f (c) _ src^i) df 
_€j* = f 8 * + aX ^ 1 i;O a _ " c(3~c a ") /ax + ax\ (53c) 

axap " \ax aV 2 "/^\ a ~ 3( c a ~i) " = 2_?A-_?ir 

\lc' + c a (3-c B ) -■— " 

where A and \x are defined by (53a) and (53"b),_and K 
is obtained by substituting the expression for f(c) from 
(44b) : ' 

■Z = K(X.n) = K(X+p) = K(c) = — - 3 (1 " c3/2) (53d) 

^ 3 V(3-c a ) y (c a -l) 3 

The numerical values jbr K are collected in table II. 

The lines- X = 'constant t and p = constant are char- 
acteristice since we had- so chosen the transformation for- 
mulae (51) . If, after the transformation, X and M> are 
plotted as rectangular coordinates (fig. 11c), it appears 
that the normal form (53c) of the hyperholi c ..equat ion has 
as characteristics, the sets of parallels to the A and 
Mj axes. For equation (52), which is also of the form 
(32), A = 0, B i= 4, C = 0, and the variables X and T 
are now a and M>. These substituted in the general equa- 
tion (38) of the characteristics, give: 
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a* an = o 

The two solutions of this differential equation are: 

A = constant 

and 

V- = constant (fig. 12) 

The solution Z of the differential equations (32) 
and (52) may be determined if, along a general curve, an 
element strip is prescri'Ded as boundary value. This curve 
may not, however, be a characteristic. But if it is made 
up of two characteristics .of different families, it is sur- 
prising that a solution of the differential equation may 
still be determined. For this purpose, the function Z 
alone is sufficient as boundary value while no elementary 
strip may be prescribed since this would be imposing too 
many conditions. 

Let the values Z = Cp(X,n 0 ) - Cp(X) and Z = \J/( X Q , n) = 

with cp(X 0 ) = ^(u D ) bo ?iven along two segments 

A 0 A 1 and A 0 A 2 of tmo characteristics (fi<. 12) . Along 

A 0 A a there is therewith also .~iven ciZ/oH, but dZ/Bu- 

is assumed not to be prescribed; similarly, along A 0 A L . 

It is to be - observed that no elementary strip is prescribed 
along A^AqAq of Z but only the values of Z itself. 

By the method of so-called "successive approximation," it 
is then possible to find a solution Z of the partial 
differential equation (52) for the entire region AiA 0 A a Aj , . 
which assumes the given values of Z along Aj.A 0 A a . 

As a first approximation, Horn (reference 10) 

Z tt = <P(* ) + ^(H) - <P(*o.U 0 ) 

for all values A and |i of the region A 1 A 0 A 2 A 3 . On 
the boundaries AqAj^ and A 0 A a Z a becomes equal to the 
prescribed values. 



*The proof will not be ;iren here. It is carried out by 
J. Horn ( reference 10) , 1913, sec. 30, pp. 164-169. For us 
it is of importance to know only that the prescribed func- 
tion Z(X,ii) satisfies tho boundary values and the hyper- 
bolic differential equation (52) . 
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We now form with the right side of equation (52) !■ 




where the integration la to be taken over the ■ doubly hatched 
rectangle. Proceeding in thi a manner, we form ■ • 

= - J J- (a ^2z± + D ^zi + cZ^) dXd^ 
*0 ^o 

Setting 

Z(X ,m.) = Z a + Zp + Zy + 

then this sun ia the required aolution and it converges, 

aa shown by Horn, in the rectangle L^k^k^L^ . 

There will now bo shorn a laat property of the charac- 
teristics - the most important for tho application to 
shooting water. At tho aamo time, in addition to the meth- 
od of solution of (32) by series development and tho method 
of successive approximation, we- shall become acquainted 
with tho method of integration of Riemann. 

We denote "by W(Z) the most general homogeneoua lin- 
ear differential expression: 

N(Z) = a Zj X + 2b z Z y + o z yy + d z^ + e Zy + f z (55) 

where the coefficients A -to F depend only on. the free 
variables X and 'Y. The general linear homogeneous differ- 
ential equation of the second order ia the equation (32) : 

N(Z) = 0- (56) 

To the expreaBion N(Z) another one M(W) i,s made to cor- 
respond, • having -the< same coef flcient.a. .A, .B.C.etc' as in 
(55 ) , where 

M(W) 2 (AW) IT +2(BW) XT +(OTr) TT -(3)W) I -(EW) T + F ft (57) 

i • ■ 

■ 

= M(W) = A W IX +2B Wjy+0 W TT +2 V z Uj+B T - i D) + 



+ 2 W T ( B Y +0 Y-i E) + W (A XT +2Biy+0yT-D I -EY+F) (57a) 
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M(W) is then denoted as the adj**fltt of -N<Z) and the equa- 
tion 

M(W) = 0 (58) 

the adjvaist differential equation of N('Z) = 0*. Z and W 
are functions of X and Y: Z = Z(X,Y), W = T7(Z,T). MCA) = 

0 has the same characteristics as N(Z) = 0, since in . 
(57a) and in (55) the coefficients of the partial deriva- 
tives of the second order are the same and since, accord- 
ing to (38), the characteristics depend on these coeffi- 
cients only. 

3y addition of the identities: 

AWZ XI - Z(AW) IX = ^ [^AWZ Z - ZCAW)^ 

BtfZyy - Z(BW) XY = ± j^BWZ z j - ± j^Z(BW) Y j 

BWZ IT - Z(B10tv = ± j^ffZyJ - ^- fz(BTT) z lj 

CT7Z 7y - Z(Gff'; TT "= J~CT7Z T - Z(CT7)y~|' 

i- 

DtfZ-y + Z (DTT) t ='^.| DST71 

oS i j 

E¥Z T + Z(EW)y = g- |~3ZffJ 
FT7Z - ZFW = O 

there is obtained the identity: 

W N(Z)-Z M(W) = £± j^AwZ z .- Z(AI7 ; z +BT7Z Y -Z(BTT) T +DZt7^ + 

L r 

+ j^BtfZ 2 - Z (B¥)j+GTTZy-Z ( Gtf)y+EZtf (59) 

Denoting for a moment the two expressions in brackets by 
E and Q, respectively, the above equation reads: 

TT N(Z) - z Y(W) = 8P/£>I + BCi/dT (59a) 

This equation we shall integrate over the region F of the 
X,T plane; Let the boundary of the region of integration, 
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to be more- definitely fixed later, "be 0 ■ (fig. 13 ) 1 

'I .. ■ . 



) ax dT 



JJ £w n(z) - z M(isr) J ax dT = yT'OPisx -i- ao/di 
(r) .. " On 

The right side may ."by integration "by parte be converted 

into a line integral. ' There is obtained: 

JJ* |w CT(Z) - Z M(W) J did! =jf (P dT - Q dl) (60) 

(!) U J (0) 

The generalized Green's theorem (60) will pow be applied 
to the normal form (52) of the hyperbolic differential 
equation. Por this purpose there is to be set in (60) 
'A = 0, B = i, 0 h 0, D = a, E = b, and I = C. In 

place of X and Y, we have >^ and \t. The expressions P 
and Q then become: 



P = i(W Z^ - Z W^)* a Z W 

Q = i(* - z ^) + b z" w 

Green's formula (60) now reads: 



(61a) 



Sf [ W N ^ Z) " Z M ^l' d * dlJL = ^ (P dn - Q d?0 (61b) 

(P) (0) 

With thiB formula we may now prove the following I 

If Z ia a function of \ and V>* Z = Z(X,iO, which 
satisfies the hyperbolic differential equation (52) and 
for which, along a' curve from Ai to B l (fig. 14) •» which 
thus ,^ in general, is not a characteristic — an elementary 
strip is given; then by these boundary values and the dif- 
ferential equation, the function Z 1b determined in the 
characteristic rectangle AiOj^OjL 1 , which oontains the 

curve AiBj with its. end points. . ' 

In order to ehow thiB we apply the formula (61b) to 
the region G and its boundary AOBA .of figure 14, where 



*Along Aj^ therefore Z ■ and the slopes &Z/b\ and dZ/Bp, 
are ?iven where naturally along Aj,Bj , the condition dZ = 
Z^ dA + Zp, d|A must be satisfied. 
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0 is an arbitrary intericnr point (X = p , \J> = q.) of the 
characteristic rectangle AiOjB^O! 1 . In integrating along 
OB, only P dp. contributes anything; Q dX does not con- 
tribute anything, since dX = 0 . Similarly, 

0 0 



J* (p dji - q dx) = - J Q ax 



A 



since along AO H = q = constant, so that dH =0. We 
thua obtain from (61"b) applied to the hatched region G 

B OA 

JJ^ N(Z)-Z ld(W) J dX d|x = j* P dn - J q, dX + J (P dp.-Q dX) 
(G) ~ 0 A B (62) 

Non from '(61a), if the first term is integrated by parts* 
B B 

J P dn = f (* v -| - i z f + aZIT) dl , = 
0 0 b 

= 4 (W Z ) B - \ (W Z) Q - J* Z (dW/3H - a ¥) dli (a) 

0 



Similarly, by integration hy parts of the first term 
0 0 

-f * ax - + .y (- i w || + 4 z ' §* -■»») dx 

A ■ A 0 

= - £(Y Z) " + £(ff Z) A +y Z(5Y/SX' - h w) dX 

A 

With expressions (a) and ("b), formula (62) becomes: 



B B B 

J \ w g dn = *cwz) I - i z |5 dn 
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(Z,)-Z M(T)J dX dn=- (W Z) 0 + £. .[(* Z)a+_(* Z >b]_ + 
0 . B 

+ Z (3W/9* - "b W)dX -^Z(3W/aii- a T)dn + 

B u * 

+^(^11^- P d|i)' (63) 

A 

We now choo se for each point. '0 which is given by the 
coordinates X = p, p- = q, a definite function W of the 
coordinates A and n: W = W(X t li). In this function.p 

and q occur as parametera, the function W(X,ii) "being 

different for each choice of the point 0(p,q). We thus 

have : 

W = W(X = W(X ,u; p,q) 
where the funotion is to have the following properties': 

1. At the point 0 itself (p,q), W is to assume 

the value 1 . 

2. The function W 1r to satisfy over the entire re- 
gion G (fig. 14) the adjunct differential equation M(W)= 
0; i.e., be a solution of 

M(W) = 0 (64) 

3a) Along the ■ straight line OB (X = p constant, 

ll variable) the function W is to assume the 

values : 

U . 

/a ( p , p ) dii 

T(p.|i) = e q (65a) 

Condition 1 is thereby satisfied since for the point X = 
Pi H = q, V(p,q) = e° = 1. Differentiating (55a) with 

respect to \i t there is obtained for the function W along 
OB the relation ■ ■ 

aW/dn - a W = 0 (66a) 

3h) similarly along the straight line AO (n=q con- 
stant ;X variable) the function is to assume 
the values : 



44 



N.A.U.A. Technical Memorandum No.- 934 



W(X,q) = e p 



(65h) 



Here, too, «lie condition UCp.q) = 1 is satisfied. Differ- 
entiating (65b) along AO with respect to ^ there, is oh- 
tained along this line the relation: 

a\t/dh - b = 0 (66h) 

The function defined by the conditions 1, 2, and 3, 
is known as Green's function W^,p;p,q) of the ■differen- 
tial equation N(Z) =0. It is determined only by the co- 
efficients of this equation. That it exists we know for 
ffi according to condition 2, is a solution of the partial 
differential equation of. the second order (M(W) = 0, for 
which the values of W alone; the two characteristics AO 
and OB are prescribed accordinq to requirements 1 and 3, 
as boundary values. It is thus poRsihle to determine W 
by the method, for example, of successive approximation. 

Substituting now in (63) U(Z) = 0, and Green's func- 
tion W, with its properties (64) and (66a, 1)), there is 
obtained : 



B 



i[(ira) A + (rz) B j+y (q d* - p dn) 



0 = - ZO + 
so that ""A 



B 



ZO 3 Z(p,q) = & |:WZ) A + (T7Z) B + (Q d?^ - P dp.) 

I A s 



:67) 



Substituting further the expressions (6la) for P and Q, 
we have : 

r 

zo = z (Pi ( i) = i ^ wz )a + (WZ)B^ + 



B 

/ (i wzx -i" zw x +hZw)a*+(- iwz^+ i ZV^-uZff) dp. = 

3 i C ) A + ( WZ ) B J + I* \ £ TT(BZ/aX cos cp - 6Z/8p sin qp) 

/ A 
B/ 

[ - £ z (Btf/bX cos cp - off/dp sin cp) 

dp 

J + z w (h cos cp - a sin cp)i as (67a) 

A/i* J 
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We here thus expressed the required solution Z at 
"point ^OCp.q.) "by -the-.given , boundary, values; itC, by a 
portion of the elementary 'strip AiBi. The c'onsi derations 

hold for every arbitrary point 0 which belongs to the 
Characteristic reotangle determined by the points Ai arid 
B x , It may be remarked further that Z Is already deter- 
mined at point 0 "by its elementary strip along AB and 

therefore that the portions AAj, and BBj, (fig.* 14) of 

■the boundary value strip A x B! have no effect on the 
value of Z at point 0. ■ " . ■ 

By means of the elementary • strip A^Bi therefore, 
.the solution Z(X,n) of the differential equation BT(Z > = 

0 is certainly determined In the largest characteristic 
rectangle which is fixed by A 1 B ] _. We wish to show, fur- 
thermore, that it i s determined only within it, and not 
outside of It. Let Q be a point without A3.O2.BiO]. 1 . 
Z is- not determined in Q since, according to formula 
(67a) ZQ depends on the elementary strip AH (fig. 14). 
The portion Bj_ H of this required elementary strip, how- 
ever, is not given. Thus the above theorem is proven. 

A special case nhich we still must examine in partic- 
ular, is that for which the curve Aj^Bj - along which an 

elementary strip of z Is s;iven - degenerates into the 
line AiOx'Bn (fig. 15), consisting of two characteristics. 
From the method of successive approximation, we know that 
Z is then determined in the region Aj.OiBj.0] 1 by the as- 
signment of the values of Z alone, along BjOj. 1 A . This 

fact will now also be derived from Rlemann's method of in- 
tegration. 

We start from the solution 

B 

Z(p.q) = * | (W Z) A + (W"Z) B + (Qd^-P dn) (67) 

U IS 

(AO x »B) ' 

Since along ±0 1 l X= constant, dk = 0, and along 0j 'B 
M- = constant, dn - 0, the integral on 'the right side 
breaks up into two-part integrals 
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B _ 0 X 1 B 

J* (Q - P dn) = J* - P dj* + J* Q 



A-O^-B A °i' 

Substituting in the above the expressions P and Q (equa- 
tions 61a), there is obtained, as before: 

Ox' A 

" J P(Ut = + y^ < 8 w ~ dZ /*v - i z aw/an + a z w)dn 

A oi' 
This time me integrate the second term by parts and obtain: 

- / p in = J(if z) 0i ,-i(Trz) A + / w(az/an + a z)an 

A Oi' (a) 

Similarly (again the second term integrated "by parts): 

J* q = -£(w z) q , - j(i z) B + y° w(az/ax + * z)ax 

Oi» 0; ' (d) 

Substituting (a) and ("b), we have, finally: 

Z(p,q) = (W Z) Qi i + J W(dZ/a|J.+a Z)dn + W(dZ/oX+D Z)dX 

Ox' Ox' (68) 

With the prescribed values of Z as boundary values 
aZ/an is also given along C^'A.- The integral from 0 X ' 

to A nay thus be evaluated without the necessity of 
giving also oZ/dX and honce p.n elementary etrip. Simi- 
larly with the Z vnluos alone , the values dZ/dX along 
0 r , 3 1 and also the second integral in (68) may be evalu- 
ated by assigning Z alone. The formula (68) thus repre- 
sents the solution Z(p,q) in the entire characteristic 
rectangle 0 1 B x Ox 1 . 
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9. Summary ?- 

Trom the differential equation of. the .velocity poten- 
tial (15) of a compressible flow and from the flow space, 
w'e were led by the Legendre contact transf ormatio'n to the 

differential equation of the poeit lon-determining poten- 
tial X (31) in the velocity plane. In aonneation with 
this partial differential equation of seaond order, we be- 
came familiar with the characteristic curves and some of 
their properties. For 'shooting water and for supersonic 
flows, these aonsfst 'of two real families of curves,' name— 
ly, epiayalolds . The Riemann method of solution showed 
that the solution of the hyperbolic partial differential 
equation by the boundary values is always determined with- 
in a complete characteristic reatenqle, namely, the smallest 
rectangle which contains all the boundary values. 



THE METHOD OF CHARACTERISTI CS 

10. Introduction 



Important contributions to the solution of the differ- 
ential equation of two-dimen b1 onal supersonic flows have 
been made by Prendtl, Meyer, Steichen, Ackeret, and Buse— 
mann . Whereas the first solution methods are purely com- 
putational,, it was pointed out by J. Ackeret that, with' 
the aid of the characteristics a graphical method may be 
developed. This has been carried out for flows without 
energy dissipation by Prandtl end Bnsemann. Tor the case 
of flaws with Impulsive discontinuities, Busemenn has de- 
veloped - on the basis of the method for nondi ssipative 
flows — a graphical method where the characteristics are 
replaced by the so-called "shock polars" (references 1 (or 
2). 7, (pp. 421-440), 14, 15, 17, 18 (pp. 499-509), end 27). 

Let the velocity of a two-dimensional supersonic flow 
or a shooting-water flow be given along a portion of a 
curve AB (fig.' 16) . Let the flow be from left to' right, 
0' a point downstream through which .pass the two Mach 
lines BO 1 and .AO 1 . . The re gion of the flow bounded "by 
the Yeah lines .OA, OB, BO', AO', we shall denote as the 
Mach quadrilateral, We shell assume that no restriction 
of the flow (vertical walls) is located in its interior:: 
that is, neither boundary nor any other object. It may be 
shown by a simple consideration' that under these assump— ■ 
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tions the flow, if prescribed along AB, determines the 
condition in the entire Mach quadrilateral AO'BOA. Out- 
side of this quadrilateral, influences from other points 
are effective. At point F, for example, another wave 
GT may arrive and pro duo e a disturbance without producinq 
a change on- AB, since. GT is a wave of the same family 
as BO' . 

Since every nondissipative flow ie also a possible 
flow in the opposite direction, the same considerations 
apply to the upstream region AOB. This statement is not 
in contradiction of the general fact that in a flow with 
the above critical velocity, the effects of disturbances 
make themselves felt only downstream. We aO not state 
that the condition at 0, for example, is caused by effects 
on A3, but rather, from the effects on AB, conclude as 
to the upstream-lying causes. 

It is to be observed that the Mach quadrilateral 
AO'BOA in general has curved sides which, as Mach lines, 
are determined with the flow itself. In the preceding 
section, from the integrals of the hyperbolic differential 
equation, we became familiar with the remarkable fact that 
boundary values act as determining factors only within re- 
stricted regions. To the characteristic quadrilateral, 
the region of solution of the differential equation, there 
corresponds in the flow the Mach put. dri lateral . The Mach 
. lines are no other than the "characteristics" of the dif- 
ferential equation of the velocity potential. The charac- 
teristics in the flow plane are 'not given, however, in 
advance as those in the hodograph, but become known simul- 
taneously with, the solution $(x,y). This is due to the 
fact that the coefficients of that partial differential 
equation (15) contain not only the free variables but also 
the. first- •derivatives of the function <D , that is, $ x 
and $y. This is-also the reason why we passed from the 

. flow space to the velocity plane (equations (31), (31a), 
and (53c)). 

IT, Physical Basis of the Method of Characteristics 

. By means of the characteristics in the velocity plane, 
it is simple to draw the field of flow of two-dimensional 
supersonic flows and also shooting water if the flow of 
approach and the side boundaries are given. With a veloc- 
ity prescribed alone; a line, the flow may be determined in 
general in the circumscribed Mach quadrilateral.' It is 
thus a question of Graphical metho'd of solution of the par- 
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tial differential equation ( : 15) or- (31). The flow is known 
"if- the velocity- ..(il,t) is known at each : point (i,y)» 
■Hence, it is not necessary to know "the Telocity potential 
4>(x,y) or- the • position-determining potential X(u,v) them- 
selves. It Is euf f ioient . only to determine Xu*Xy 'and 

*X"*y (Compare formulas (29): Xu B x, Xv = Y and 

* x = tU <by m V. ) 

1 The graphical method Is based on'the simultaneous. con- 
struction of .the flow In the velocity field (u t v) and In 
the field of flow (x,y). , . 

Letus consider first a parallel-flow assumed to be 
bounded on one side. At the position S, the wall re- 
ceives, a small deflection 8 (fig. 17). In the case of 
supersonic flow and shooting water, this leads to a pres- 
sure increase.* 

If the wall has a convex corner, a flow arises with 
diverging cross section. In the case of shooting water, 
this leads to a level drop and acceleration. * 

Since in the boundary of the frictionless flow of 
figure 17, no finite length occurs as reference length, 
all streamlines must "be similar with respect to the corner 

8. Water depth and velocity in magnitude and direction 

therefore have constant values along-each stream through 
the corner. 

The flow of figure 17a for large' deflection angles, is 
described in Part II of this report (T.M. No. 935), under 
Shock Polar Diagram, page 1. This flow is nonstat ionary . 
The discontinuities of the different streamlines are equal 
and all lie on a straight stream ST passing through the 
corner. For extremely smnll deflections, the corner leads 
to only a small dlsturbnnce in the flow. Since small dis- 
turbances have the Mach linos as the wave front, the dis- 
turbance line ST is a Mach line. It forms with the 

*The following considerations hold for 'water and gas flows. 
Since, however, for tho analogous concepts different terms 
are applied in hydrodynamics and gas dynamics, both would 
always have to "be carried along In this work. This di.ffi» 
culty has.heen avoided as far as possible hy using the 
terms from hydrodynamics. Where terms from gas dynamics, 
nevertheless, occur the corresponding terms are: Expan- 
sion = level drop; compression = level -.rise; Impulse = 
jump; oxpansion .wave = depression wave, etc,.. 
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parallel flow an angle 'a where sin a = a/c = 
For somewhat larger deflections the discontinuity lies on 
a stream ST, Whose direction lies between the directions 
of the two Maoh lines of flow I "before the deflection , 

and flow II after the deflection. 

The flow corresponding to figure 17b for large de- 
flections and hence, strong acceleration, Is treated more 
In detail In section 21, Part II of this report (T,M. No. 
935), under Level Drop about a Corner. In contrast to 
level rest, the drop is continuous. It "begins again on 
account of the similarity for all streamlines on a stream 
ST'. This Is a Mach line of flow I before the level drop. 
The deflection for all streamlines ends on a stream ST" , 
a Mach line of flow II. For small deflections, It may be 
assumed as a first approximation also for the level drop 
that It Is concentrated on a mean stream ST. An impor- 
tant simplification Is thus obtained for the qraphlcal 
method . 

Both the small level drop (In the gas expansion) and 
small level rise (compression) have the follomlnq in com- 
mon: The velocity receives along a disturbance line a 
change in magnitude and direction. The direction of the 

disturbance line is given as the mean direction of the 
two Mach lines of the conditions before and after the 
change.* In traversing this line, there Is also a change 
In the pressure. The pressure drop or gradient - that is, 
the Increase In pressure per unit length in the direction 
of the most rapid change - is thus normal to the mean Mach 
line. According to Newton's law, the acceleration and 
hence also the vector change In th'o velocity, has the di- 
rection of the force. We thus have the result: The ve- 
locity vector Cj before the deflection (rise and drop) 
receives as a result of tho deflection, a vector Increment 
Ac which is normal to the Mach line. Since the deflec- 
tion angle Is also known, Ac Is determined (fig. 18). 

The graphical method consists In building up the en- 
tire field of flow out of small individual Mach quadrilat- 
erals, in each of which the velocity is constant and de- 
flections occur from ono quadrilateral to the 'other. 



*Wherever necessary for clearness In what follows, a dis- 
tinction will be made betweon disturbance line and Mach 
line. The disturbance lines are those along which the dis- 
continuities arise.. Disturbance lines of Infinitely small 
intensity are - Mach Tines. Both pass over into one another 
In steady flow. 
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12. Mach Number and Angle. ■. ■■ 

It IB- important that the Maoh'hum'ber " If : and : the angle 

a (sin a = l/lf) are given by the magnitude of the flow ve- 

-locity alone, since sin a = J~kh./o and. according to the 
energy equation', the water depth h depends uniquely on 

the flow velocity (equation (9)) . We thus have: 



Biu 3 a= «h/c 8 = (gh G - £ c a )/c 1 

m 

(42) 



o 

Dividing numerator and denominator of the right side by 

a* a = 2gh Q /3 

we obtain in the notation of nondimensional velocities 

C = c/a* : 

l/M 8 = Bin a a= ("I -| c a )/c a (69) 

For the graphical method, there ia applied the graphical 
representation of equation (59) (fig. 19 ) ■ a being plot- 
ted as arc, _ and >c,'as radius vector. In rectangular co- 
ordinates, v = c sin a, 

— a — a , a s 1 — 8 

v = c ain a, = -r - ■= a 

and 

_a m * a \ 3 -a 3 

u =c (l - sin a; = 2 0 "* 2 

Eliminating "c from those two equations, there Is obtained 
the curve in rectangular coordinates 

v 3 =1 (70) 

This is an ellipse with major nnd minor semiaxes n/"3 and 

1 (fig. 19) . For an ideal 3.1.8, it is an ellipse with the 

semiazea J (k + l)/(k -. l) and 1.. 



13. Oh.iracteristics 

If any nondimensional velocity. Cj .is given at point 

P of the flow .plane, the direction of the Mach line at 
the point considered is obtained in the following manner: 
"5j is drawn, in the- velocity plane (Jig. 20). The .ellipse 
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is now rotated about 0 until the extremity of Oj lies 
on it (two possible oases) . . Then, according to figure 19, 
the principal axis of the ellipse so rotated gives the di- 
rection of the Maoh lines in the flow and according to 
figure 18, the minor axis of tho ellipse gives the direc- 
tion of the velocity increment Ac. Pour typos of increase 
are possible, depending on whether the Mach line is a dis- 
turbance line of the first or second family, and whether 
the disturbance is a drop or p. rise. In the example shown 
(fig. 20) no disturbance line of the first family passes 
through the point P, whereas that of the second family re- 
sults in a deflection, namely, a level drop. The velocity 
'increment, denoted by a heavy arrow, thus, is the one that 
comes under consideration for this example. If the dis- 
turbance lines of both the first and second families pass 
through the point P, the apparent difficulty is removed 
by considering a neighboring streamline. For the latter, 
the velocity receives tao ohanges, one following shortly 
after the other, each of which is uniquely determined. 

At each point of the velocity plane there are thus 
two directions of the velocity Increment. These two di- 
rections are given by the minor axis of the ellipse (figa 
21.* There is thus obtained in the circular ring area, 
between R = ^3, and r = 1, a direction field which de- 
termines two families of curves. In figure 21, two repre- 
sentatives of theso two fnmiliofl are drawn. By the fol- 
lowing simple consideration, Buooinann shows that we have 
here the case of the previously found epicycloids. 

The direction field is obtained by drawing the small 
segments a, b, c, d, ... in the direction of the minor 
axis of the ellipse (0, , l), then rotating the ellipse 
some-hat, and again drawing the linns. We may now consid- 
er a, o, c, . . . as lying, instead of on the ellipse, on 
the fixed points of the circle chords AiAq , Bj.Ba , Cj Cg , ■■■ 

There is then obtained the same direction field as_ before 
if these chords are rotated in the circle (0, a/iT) and a, 
"b, c, . . . drawn each time. If all these chords with their 
points a, "b, c, . . . are now arbitrarily drawn in the cir- 
cle (0, J~2>) (fig. 22), the small segments a, t>, c, . . . 
are still in the direction of the required direction field. 
By suitable rotation of the chord diagram (fig. 21 ), we 
pass a family of ohords through an arbitrarily chosen point 
A i, the chord diagram "being rotated so that B L , 0j_ , H 1 , . . . 

""Figs. 21, 22, and 23 correspond to figs. 40., 41, and 42 
of Busemann, 1931, p. 422 (reference 7). 
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lie successl vely on Aj. end the segments 8, b, •-■ ■ 

Being • drawn. ■ Ehe latter- *rl 10.- st ill. he segments, in the 41- 
reotion field (fig. 23); -The complete field will '"be OB- 

tained by rotating thie diagram about 0; 'for example,' 

Ax toward Ax 1 , and then again drawing the small segment a 

a, D, c, . . . 

Now the points a, 1), a, ... . divide the chords Aj.Aa • 

B i s a • °i ... (fig; 21) in the same ratio; the ellipse 

as effine figure of the circle having this property: The 
points a, "b, c, ... in figure 23, thus lie on a circle. 
The directions a, b, c, . . . are normal, respectively to 
Ab, Ac, . . . 

If the circle with diameter AAi is rolled on the' 

circle about 0 with the radius 1, each of its points de- 
scribes en epicycloid. The rolling circle at the instant 
represented, rotates about the point A. All of its points 
thus also move on normals to the lines joining the corre- 
sponding points with A, the direction field of the set 
of epicycloids Being identical with that of the required 
curves of the possible velocity Increment Ac. These carves 
are thus the epicycloida described above (figs. 21 and 9)» 

We have mentioned the same epicycloida before. They 
are the characteristics of the partial differential equa- 
tion of the flow. We now see the physical interpretation 
of the characteristics: Curing the passing of 8 smell 
disturbance wave the flow velocity changes along the cor- 
respdnding characteristic. 



14. Graphical Construction of the Tlow 

The field of flow and the hodograph are drawn simul- 
taneously - .in the hodoyraph, the velocities and their 
changes; in the field of flow, the streamlines* The flow 
is always assumed from left to right.' TTe may then speak 

of en upper or a lower boundary. All disturBanoe lines 
that start from the upper boundary will Be. denoted as the 
upper system of waves, nnd'all thoee from the lower bound- 
ary', 'the lower system. 1 

a.) Flow bounded on one aide.- The simplest supersonic 
flow is that bounded on only one side .as given by the bound- 
ary conditions of figure 24. Let the' approach be parallel 
end have the Mach number II = 1.5. ' As a first step the 
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continuously curved mall is replaced by small straight se%- 
menta with angle Increment e of, . for example, 2* , In some 
Cases it may be Of advantage to make the angle increment s 
of various amounts. 

To the flow of approach (parallel flow) , there corre- 
sponds, in the velocity plane, a single point Pi given 
by the direction of Oj . and the magnitude Ci*. Pi is 

also obtained as the point in the hodograp.h (fig. 24c) at 
which the normal to the characteristic forma with the ve- 
locity, the Mach angle cti • At the flow receives a 
first discontinuity, a level drop which leads to a deflec- 
tion by the angle 6. This deflection is of equal magni- 
tude for all streamlines and lies for tho entire flow along 
the disturbance lineSj^, whose direction wo shall learn 
from the hodogrnph. In the latter the velocity c a aftor 

the first discontinuity is given by the point Pa,, whose 
radius vector forms the angle 8 with that of P x , and 
which lies on the characteristic through P L , correspond- 
ing, fo-r c 1 c 2 , to a drop; that is, an increase in 

velocity. We thus obtain P a and c a . The disturbance 
line Sj_ Tj_ in the flow is, as ire know, a mean Mach line 
between the states Pi and Pg . This direction is now 
given simply as the normal to the characteristic between 
Bl and Pg in the velocity plane. In the entire region 
2, the flow is again a parallel flow with the Velocity Cg 
up to the disturbance line S a T a . This line and the stat8 
after this second di at urbane a, i s determined similarly as 
for Si Ti , only now the initial velocity is given in the 
hodograph by P a . The velocity after the disturbance is 
again the velocity OP, deflected by 6. The direction 
of the disturbance line S a T 2 is the direction of the nor- 
mal to the characteristic between P 3 and P 3 , etc. 

Withthe above construction, the first disturbance 
thus lies along S 1 1E ll the last along Sn^Tn-i. Actually 

the "beginning and end of the disturbances lie along the dot- 
ted lines S 0 T 0 and ST, which have the directions of the 
normals to the characteristic in Pj_ and ^ . It is only 

. m 

*Prom equation (69), we have: c a = 3 M a /(M a + 2) 

For gases: c a = (k +l)li a /[(k- l) M a + 2] 
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beoause we must draw the flow disoontinuously in finite 
■*bep«- that- the aotual-....atart : of the disturbance and the 
f ir.at. diaturbance do not accurat ely coincide BP .de'Creaa- 
ing- the steps, 'the accuracy may be raised. 

Figure 25 shows a flow drawn in this manner with M = 
1.5, and for rater (t = ,2), the deflection increments 
"being 2°. From this simple example, an important property 
of shooting water "bounded on one aide (supersonic flow) 
nay be recognised, namely, that as' long as no large discon- 
tinuous pressure riaea (impulses) occur, all the pb-ints 
giving the state in the hodograph lie on a single charac- 
teristic; i.e., for such a flow the magnitude of the ve- 
locity depends uniquely on it a direction and viceveraa. 

A limiting case of the example considered is the level 
drop about a corner (fig. 26a-c) (references 14 and 17). 
This flow is a parallel f low with a Mach number equal to 
or greater than one. The one-aided rectilinear boundary 
ends at S. On thd lower side of the boundary the water 
depth (pressure in the gas) is zero or at least smaller 
than in the parallel flow of appr oe c i. The same results 
hold as for the - flow nf figure 24 except that now the 
lines SiTi.SgTa,. . . all pass through the point S . ■ The 

velocity varies along a streamline in 'such a manner that 
its end point travels on n, characteristic in the velocity 
1 plane (fig. 26c). The constant velocity along a stream 
SP has its end point P' at that iDosition of the corre- 
sponding characteristic where the normal to the character- 
istic is parallel to SP. 

b ) I nterior of a flow bounded on two s ides *- Let the 
velocity o x be given in the interior of a flow in a cer- 
tain region 1 (fig. 27). Let this region be 'bounded on 
the right aide by an upper (b), and a lower, disturbance 
line (a). -The streamlines rt and 0 , which-may 'also be 
considered : as. walls, are correspondingly assumed to have 
small defleotiona at A and B. ■■ The deflections' 8 a and 
Bp are given. The point P, in" the hodograph is the im- 
age point of the region 1 of the flow (fig. 27b). In 
crossing the disturbance wave a from region 1 to region 
2 (drop, since def lection i a toward outside) the velocity 
Ci receives a change such ' that the velocity c a Ilea on 
the characteri stic corresponding to the lower disturbance 
wave system and forma with "ci the angle 8,. This gives 
the point 'Pa. in the hodograph as in a flow bounded on 
one side and hence also the direction of a as normal, to 
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Pi Pa . The same la true in crossing the disturbance wave 
b. To this corresponds in- the velocity diagram a travel- 
ing along the characteristic of the upper system from Pi 
toward P3 (6p is given) . At a position X the two 
disturbance waves meet and their effects will "cross." 
From the point X a disturbance wave of the lower set a' 
starts out and one from the upper set "b 1 . Crossing 'a 1 
In the flow means in the hodoqraph, as in a flow bounded 
on 'one side, a change in the velocity from P 3 toward Q 4 
(fig. 27"b) where Q 4 for the present, is unknown. Simi- 
larly the velocity on crossing h' receives a change_ from 
to .S4 where S4 similarly is for the present, Tin- 
..known. Bow a first condition for Q 4 and S 4 is that 
the velocity in the region 4q of the flow on passing from 
from 1 — »> 3 — * 4 , should have the same direction as the ve- 
locity in region 4s on passing 1 — > 2. — *■ 4 . This means 
in the velocity diagram that the points Q 4 and S 4 must 
lie on a straight stream through 0 : 0S 4 || OQ^ . There 
is, furthermore, to be satisfied, the condition that the 
water depth (pressure in the gas 5 in the region 4q must be 
the same as in 4s. As long as the flow is free from im- 
pulse, the water depth is uniquely determined by the ve- 
locity. The requirement that the depth ahould be the same 
in 4q and 4s, means therefore that the velocity 0S 4 must 

have the same magnitude as 4 : 0J> 4 = (K^ . Both condi- 

tions are simultaneously satisfied if S 4 and Q4 coin- 
cide at the point of intersection P 4 . The entire region 
4 of the flow is thus in the velocity diaqram given by the 
point P4. We may now draw a' and "b 1 . They start from 
.X in -the direction of ths normals to P 3 P 4 and Ps P 41 
respectively . 

Figure 28 shows the intercrossing of two streamlines 
where now one disturbance is a level rise, the other a 
drop. The picture would be quite similar if the two dis- 
turbances were level rises. 

We shall now follow a disturbance line in the interi- 
or of a flow in the case where it, encounters several dis- 
turbance lines of the other family (fig. 29) . 'The direc- 
tions of a, h, a', and T> 1 and the points. Pj. , Ps. P a , 
and P 4 are assumed to be determined by the method given. 
Then for the regions 3, 4, 5, and 6, we again have P 4 
and Ps lying on the characteristics through- P 3 . The po— 



W.A.C.A. Technical' Memorandum Nq. ^934 



57 



sitlon of P B 1b determined by the deflection -835 and 
E* " is ■fdxed by-the charact eristics .^P*, and P 3 P 4 . 
There Is now obtained also P 6 and hence the velocity' 
0P a in region 6, Ps being "the point of Intersection of . 
the two characteriBti.es P R P 6 H and P 4 P 6 . Similarly, there _ 
is 'finally obtained P a .. The individual portions of the 
disturbance wave aa 1 a" a 1 " are In the directions 'of the 
normals at the centers of the 'portions of the charact eri s-* 
tics' P x Pg, P 3 P 4 - , P B P e , P ? Ps , respectively. 

We thus' find the result, namely, that the extremities 
'of all possible velocity vectors before crossing the di"S- 
•turbance wave aa'a" . ... the points Pj, ,P 3 |P F , ■ '• . . all 

lying on a fired characteristic through P x ■ Similarly, 
all extremities of the velocities after crossing thedis— 
turbanoe wave a - that is, the points P a t P 4 ■ P 6 • . m . lie 

on the characteristic through Ps. Crossing the 'disturb- 
ance wave aa' a* 1 a" 1 at any position in the direction of 
the flow, has the result with respect to the velocity, 
that there is ■& transition from the characteristic 1 to 
the characteristic 2 (both of 'the same family) each time 
along a characteristic of the other family. These changes 
are the heavily drawn portions of figure 29b t Since the two 
families of characteristics lie symmetrically: 

<* p 7 op 8 = * p 5 cp 6 = <$ p 3 0? 4 = <J P, 0P a = o ia 

i.e., 

81a = 634 = 6 56 = 6 78 = . . . . 



In figure 30, let the curves denoted by K be circles 

about 0.. We then have: 

a) 4 AOC = £ EOF, because each characteristic of the 

same family arises from the other by rotation 
about 0 . 

b) <fc JLOB = 4 BOE = 1/2 ■£ AOE, because A.B is sym- 

metrical to -B3 with axis of symmetry BO. 

c) £ COD = £ DOF = 1/2 £ GQP, similar to b) . 

d) ' <3t COE = ^ COE. 

Equation d) subtracted from a) «;ives 
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4aoc-± COE = £ EOF - ■£ COE 

i.e., A033 = £ COF, and hence it follows from h ) and c) 
4- ?0E .= ■$". DO! i as was to be proved. 

Tfe thus obtain the moat important result: On cross- 
ing a disturbance wave the velocity undergoes a change in 
magnitude and direction. The change in the velocity di- 
rection is the same at all points of the entire disturb— 
ance wave Independent of the direction of the velocity be- 
fore the arrival of the di sturhance wave and regardless of 
whether or not the wave was crossed by disturbances of the 
other family. This is true on the assumption of flow free 
from impulse. In section 4 we consider flows with Impulse 
for which the velocity is not a unique function of the wa- 
ter -depth. 'There it will be found that the deflection 
angle caused by a disturbance wave may vary along the wave. 

c) Fixed wall wi th 8 flow bounded o -n -two- si des » - In 
figure 31, let SAC be the upper boundary of a flow. Let 

no disturbance wave from the opposite wall meet the corner 
S of the wall at first. From the latter, 8 wave s 

starts out which is identical with that of a disturbance 
starting from a flow bounded on one side. 

We shall now consider the effect of a disturbance wave 
a which encounters the straight wall SC at point A. 
In region 1, let the velocity be given by the hodograph 
point Pj_ (fig. 31h). On crossing the disturbance wave 
a from region 1 to region 2, the velocity receives a de- 
flection 8, given hy the lower wall. Pa lying on the 
characteristic is thereby determined and also the disturb- 
ance line a. Since at each point of 'a flow there are two 
possible disturbance waves, there can start out from A 
only 8 wave of the uppor family ("b). The line "b and the 
velocity in region 3 are determined from the condition 
that first the velocities c x in region 1, and c 3 in re- 
gion 3, must be parallel, since it was assumed that the wall 
had no discontinuity at A. In the hodograph this means 
that P 3 must lie on the straight OPx ■ Secondly, "b is 
a disturbance I i ne from the family other than that of a, 
so that P 3 lies on the characteristic P a P 3 , which pass- 
es through P a . By both of these conditions P 3 , the ve- 
locity c 3 and alBO the disturbance line "b are deter- 
mined. 

The angle of deflection which the velocity undergoes 
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on cro a sine; the reflected wave la equal and opposite to 
■the^-angle. of -deXle-Ctlon by, the incident disturbance line. 
If the incident disturbance i A s' a level rise, then there- 
fleeted disturbance is also a rlBe (fig. 31b). If the 
disturbance line is a drop,_ then the reflected line' is also 
a level-drop disturbance (31c).. 

In case the disturbance line' a strikes the wall at 
the position' S where, the pall has a discontinuity, no 
new difficulty arises. It is then only necessary to imag- 
ine that the reflected disturbance line b and the newly 
generated dleturbance line "B follow shortly upon one an- 
other. If h and ' s are both level-drop waves, each 
must be. drawn separately;, if both are level-rise waves, 
then they are drawn together, as a single disturbance start- 
in? from S, on the crossing of which the velocity under- 
goes a deflection equal to the sum of the deflection8 due 
to s and b. If, however, one of the disturbance lines 
is a rise, and the other a drop, then .only a single dis- 
turbance line starting from . S is drawn* along which the 
deflection angle for the velocity is equal to the differ- 
ence between the deflection anqles for s and "b and, 
depending on. the Intensities of s and' b , may be a rise 
or a drop line . * 

In the third case, whore the deflection angles for s 
and b are opposlto, it mny'also happen that they have the 
same magnitude. In that case no disturbance at all start 8 
out from that point. This is the case if the wall itself 
has the same deflection angle as that of the approaching 
disturbance wave. This fact is made 'use of where it is de- 
sired to produce a parallel flow. In the latter no dis- 
turbance waves occur. This .condition is obtained by giving 
the walls in . succession di scontinui ties such that one dis- 
turbance wave fs "swallowed" when the other strikes it. 

d) Free .let .— If a disturbance line strikes a free jet, 
another type of reflection occurs since the water depth 
must have a fixed value (fig. 32). Let the point P x in 
the velocity diagram .correspond to region 1 ahead of -the 
disturbance wave. The point P a which gives the 'velocity 

*Por the third-case it is clear that only a single disturb- 
ance line starting from S is drawn because the sum of the 

two disturbances is smaller than that of either Individual 
case. For the first case two, and for the second base only 
one, disturbance line is drawn in order to approach the 
true condition for which drops are spread out in the form 
of a fan (drop about an edge; while rises are concentrated 
(impulse) . 
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0P a of region 2, Ilea on'the characteristic through Pj_ 

belonging to the lower'' family of disturbance lines and de- 
termined by the deflection angle 8, . Since at each point 
two disturbance waves, at most, pass through, there oan 
etart'out at point A of the flow where the line a 
strikes the free jet, at most, another disturbance line b 
of the other family (b) . ' The disturbance "b must be such 
that the water depth is the same in regions 1 and 3. This 
means for flow without energy dissipation that the hodo— 
graph point P 3 corresponding to region 3, must lie on a 

circle through Pj_ about 0: OP, = OP,. Since, moreover, 
P 3 " lies on the characteristic through P a belonging to 
the upper di sturbanee line, family P 3 is uniquely deter- 
mined and hence, also b. On account of the symmetry of 
the two families of characteristics £ P x 0P a = P 2 0P 3 . 

A level-drop wavo is reflected on a freo jot as a level- 
rise wave, and conversely. It is important to ohserve 
that the velocity deflection on crossing tho reflected save 
is as large as that on crossing the incident. Here again 
we find that disturbance waves - whether they are crossed 
by. others or reflected - produce at all points equally 
large deflection angles of the local velocities. 

15. Application: Laval Nozzle 

Let a Laval nozzle be drawn for mater (k = 2) in 
which the flow is parallel at the minimum cross section 
(M = l) and which is to produce at its exit a parallel 
flow of Mach number M = 2. 

Aside from flows with hydraulic jumps (shocks), all 
the phenomena have been discussed in detail in the previ- 
ous sections. There are no difficulties in drawing up the 
flow with the aid of tho basic elements described above. 
Instead of drawing Mach lines, however, as normals to the 
characteristics, the accuracy is considerably Improved by 
using the ellipse construction described in sections 12 
and 13. The normal to the characteristic is then obtained 
as the direction of the major axis of the ellipse without 
requiring eithor the tangent or the. normal of the charac- 
teristic itsolf (figs. 20 and 33).. 

A convenient arrangement for the drawing is shown on 
figure 34. A strip B is glued on the transparent paper 
A with the ellipse E, the edge of the strip being paral- 
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lei to the minor ' axi s. of ■ the . ellipse and' 'rotatable about 
.a nee,dle at point ^ 0 in. the origin- of the velocity' plane. 

The direction of the major- axis f s drawn with the triangle 

1 as disturbance wave in the flow. . . . 

The Laval nozzle investigated has as its boundary at 

the approach side of the flow, a cubical parabola .with 
a ehort connecting straight piece QR, in order that at 
the minimum cross section the flow, for the ahoot ing-water 
region to be drawn, ehould be parallel.. Therewill then 

be no disturbance waves in ft. To the straight portion_ 
there is connected a circular arc -RS . The shape of this 

'portion can he' chosen at will and the first disturbance 
waves start out from It. The shr\pe of ST is determined 
"by that assumed' for RS since the former must be such 
that, starting from the channel exit, there are no disturb- 
ance waves in the flop. 

If the approach flow is parallel, the construction of 

the flow begins with the first disturbance line from RS , 
the line being that of a flow bounded on one aide. The 
construction is then followed As discussed in the preced- 
ing paragraphs. 

Since we Are constantly pas sin? from the velocity di- 
agram to the flow diagram And in order that corresponding 
points may be rocognizod as such. It is necessary to intro- 
duce a suitable notstlon. For this purpose the curvilin- 
ear coordinates A And \x are convenient (equations 
(53a) and (53b)) . The numbering is shown in figure 34. 
The number beside each characteristic of the upper family 
3ives the Angle in degrees At which It starts on the unit 
circle, and similarly, for the coordinates of the charac- 
teristics of the lower family. In order that the two fami- 
lies of characteristics may not be confused, the coordi- 
nates of the upper family are preceded 'by a zero.* The co- 
ordinates A And [J. of the velocity plane Are written in 
the corresponding field of flow". The numbers thue written 
have the .property (equations (53a, And b) ) that (A - 

2 = Cp; that is, their hn.lf difference gives the angle of 
the flow with respect to the horizontal. Their hnlf sum 
(A + Ji>)/2 is A number on which the magnitude of 'the non- 
dimensional velocity and hence also the _ water-dejoth ratio 
k/h G uniquely depends, since \ + U i-8 constant' on' dir- 

*To the curvilinear coordinates ^ = 0, |* = 00, for exam- 
ple, correspond the polar coordinates -c = 1, cp = 0 . 
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cles about 0. With a definite value (A + n)/2 1b asso- 
ciated' th-e same water-depth ratio o-, JdJ* (gas temperature 

ratio T/T 0 , hence pressure ratio, p/p Q ) , which corre- 
sponds to the level drop about a corner starting from 
M = 1 (fig. 26"b) and deflected from the direction of the 
approach flow by the angle u) = ("h + kO/2. Corresponding' 
values h/h 0 , p/p o i M, O, and w = (A + n)/2 are col- 
lected in tables I and II. 

In general, the difference of the two coordinate num- 
bers is not required since the direction of the stream- 
lines in 'each field may be taken directly from the. veloc- 
ity diagram. The streamlines may also be simply and rap- 
idly drawn with the arrangement shown in figure 34, it "be- 
ing only necessary to pass the major axis of the ellipse 
through the holograph point given by the coordinate num- 
bers, the triangle then giving the velocity direction in 
the corresponding field. 

The aum of the two coordinates, however, is required 
If it is desired to draw the lines of constant water depth 
in the flow. These lines may also be drawn without know- 
ing the coordinate sum if equal deflectlone are chosen 
for all disturbance lines, namely, as diagonals of the 
Mach quadrilaterals. 

In all problems in which a parallel flow is given as 
initial flow, we begin, according to the characteristic 
method, with the first disturbance lines starting from the 
boundary . 

Under suitable assumptions, there may also be pre- 
scribed as an initial element, the velocity distribution 
along a line. The latter must not, however, at any point 
touch a Yach'line. It must thus be a line which in 'itself 
is not a Mach line and which does not intersect the same 
Mach line twice. Streamlines and their orthogonal trajec- 
tories certainly are such lines. The flow may then be 
computed by the characteristics method in the entire Mach 
quadrilateral described about this line. Thia Mach quad- 
rilateral is only determined on drawing the flow. If the 
velocity along a line is prescribed as initial element, a 
further condition is that the position of this line with 
respect to a side boundary is such that no flow restriction 
falls within the Mach quadrilateral described about the 
line except when the latter has the form of a streamline. 
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For the graphical determination of euah flows the 
line must first be broken up into suitable segments on 
which the velocity is constant- in- direction and magnitude, 
These pieces are then separated by disturbance waves and, 
starting from these, the flow may be determined with the 
Mach. quadrilateral. ■ 



List of Most Frequently Occurring Symbols 

Si acceleration of gravity.- 

$, gas constant. 

v, kinematic visoosi'ty. 

P, density. 

p, pressure. 

T, absolute temperature. 

i i heat oont ent . 

Cp, specific heat at constant pre a sure-. -- 

c v , specific heat at constant volume. 

k=Cp/c T t adiabatic exponent. 

velocity potential. 

X, posi'tioning-determining potential. 

ii7i<i rectangular coordinates in the flow space. 

r,A, polar coordinates In the flow plane (x,y). 

?Vil*i curvilinear coordinates in the velocity plane, 
charncterist lc coordinates. 

general variables. 

components of. the velocity in the x, y, and z 
directions. 

e,q>, polar, coordinates in the velocity diagram (two- 
dimensional flow) , 
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c max» maximum velocity. 

C, velocity increment. 

a, in gas: velocity of sound. 

in water: propagation wave velocity -J gh. 

a*, critical velocity. 

u, r, c, . . . , nondimensional velocities (reference velocity 
a*; in hydraulic jump a x * the critical 
velocity before the jump) . 

i£ = c/a, Mach number. 

a=( sin" 1 ) (a/c ) , Mach angle. 

h, water depth. 

h Q , total head (water depth for c = 0). 

h 0 ' ,h 0 " , total heads after hydraulic jumps. 

p Q , T Q , i Q ,li o , subscript 0: stagnation state. 

T*,h*,..,, asterisk *: critical state. ' 

Tlj , c x .hj .M^ , subscript 1: before hydraulic jump. 

u 2 , c a .i^ ,M a , subscript 2: after hydraulic jump. 

U 2?* velocity after right hydraulic jump. 

A(Z,T),B,C, coefficients of linear partial differential 

equation of second order. 

fti^iC, coefficients of the differential equation in 
normal form. 

K, coefficient of the differential equation of 
the flow In normal form. 

6, small deflection nn^le. 

UJ , deflection anejle of the flow without dissipa- 
tion (sec. 21, Part II, T.M. No. 935) . 

P, deflection an^le for hydraulic jump (fi^s. 

?S7 and 38, Part II, T.li. No. 935) . 

f , an?le of the hydraulic jump wave front (fi^s. 
37 and 38, Part II, T.M, No. 935) . 
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*See reference 7, pp. 426-7. j or values of K, see refer- 
ence 1 (or 2 ) ,p. 317. 
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TABLE II 
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tig: 1,2,3,4,5.6 





Figure 1.- Mach rays. 
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Figure 3.- Notation for energy equation. 




Figure 5.-f-«urfaca 
■trip. 
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Figure 2.- MftCh liner, . 

double family. 




Figure 4.- Skatch for 

derivation 

of continuity eguation. 



Figure 6,- Contact 

transformation 

for QU9 Independent 
variable . 
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Flga. 7,8,9,10,11,12,13 




Figure 7.- Element transformation for two 
independent variables. 



Figure 8.- Polar 

coordi- 
nates in the 

velocity diagram. 





Figure 12.- Characteristics 

of the normal 
form. Method of successive 
approximation . 



Figure 9.- Characteristics of the flow 
differential equation. 





Figure 13.- General 

region 
of integration. 



Figure 10.- Construction of the 
characteristics : 



Figure 11.- The 

vari- 
ous coordinates. 



(a) Flow plane, (b) Velocity diagram, (c) Characteristic coordinates 
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Figs. 14,15,16,17,18,19 
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Mach lines of 
lower family 
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Figure 14.- Region of 

integration for 
the normal form of the 
hyperbolic equation and 
characteristic quadrilateral, 




Figure 16.- 



Mach lines of 

upper family 

Mach quadrilateral. 




Pigure 15. e. Notation for 

application of 
formula (67) if the 
boundary value8 2 are 
given along two character; 
istics. — 
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Sida wall 5 



(a) Rise (compression) 




(b) Sink (rarefaction) Sicfe 
Mach iines, **U 

Figure 17.- Small deflection 

of a parallel flow, 




Figure 18.- The change in 

velocity on' 
crossing a Mach line. 



Figure 19.- Relation between the flow velocity 
c and the Mach angle am 
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(a) Tlow: 
p <fcV 0 Streamlines; 

Mach linat. 

(b) Velocity diagram. 

Figure 20, Imploynent of 
the hodograph 
5 for the determination of 

the Mach line in the flow. 
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Figure 24.- Flow bounded 
on one ride. 



Figure 21.- Field of directions of 
the velocity change on 
crossing s> disturbance line. 



Figure 22 and 23.- Proof that the direction field (fig. 21) 

belongs to two families of epicycloids. 
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(a) Flow plane, (b) Velocity diagram. 

Figure 27.- Interior point of a flow bounded on 
two sides '( the deflection angles 8 
which are of the order of magnitude of 1 degree 
are in thir and the following figured drawn 
exaggerated for clearness) . 



(a) Flow plane. 

(b) Velocity diagram. 

Figure 28.- Interior 

point of a <t 
flow bounded on two 
sides. 
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Fig*. 30,31,32,33,34 
Fixed wall „ 





a) Flow plane, 
(b) Velocity 

diagram for a . 
« level raiting 
(condensation) 

wave . 

c) Velocity 
diagram for a » « 
level lowering 
wave. 



Figure 31.- Disturbance wave striking a wall . 




Figure 32.- Disturbance lines striking a 

free jet boundary, 
(a) Flow plane, (b) Velocity diagram. 



Char. of loirer 
family 



Figure 33.- Sketch showing method of 

determination of the 
direction, of the disturbance wave 
by means of the ellipse. 




\_ jU = const* 



Figure 34.- Drawing of the 

flow. Increments 
of 6 degrees. 



